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SYNOPSIS 


SOBIE ASPECTS OP QUANTUM THEORY OP CHANHEIING 
/inand PrakasH Pathak 
Ph.D. 

Department of Pbysics 
Indian Institute of Technology Kanpur 
March 1971 

In the present work, the c[uant-um mechanical approach to the 
theory of channeling and hlo eking of particles in crystals has 
been investigated using Debye model for lattice vibrations and 
one phonon approximation for inelastic processes taking place in 
the crystal when the particles propagate through it. The channel- 
ing and emission of charged particles in crystals and anomalous 
penetration of neutrons have been considered to illustrate the 
various aspects of the quantum theory of channeling. The rate of 
energy loss of charged particles due to conduction electrons is 
important under channeling conditions and has been calculated 
using a dielectric formulation. 

After a review of the classical theory of Lindhard in the 
Chapter I, the need for a quantum theory has been briefly dis- 
cussed. In the Chapter IT., a general formulation for channeling 
of particles into crystals has been developed and an expression 
for the renormalization of the initial state wave function due 
to inelastic processes has been obtained using a general inter- 
action potential. The approximation involved here is that the 
ciystal is initially in a very low lying state and only the states 
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close to the initial state of the crystal are excited hy the inci- 
dent particle under channeling conditions. The general formula 
has heen simplified using first Born approximation and Dehye model 
for latticu vibrations to calculate the imaginary part of the renor 
malization matrix which governs the attenuation of the particles. 

In Chapter III, different potential models have heen employe 
to illustrate the details of the calculation. Using a potential 
for neutrons which goes over to Fermi pseudo-potential as its width 
vanishes, it has heen shown that for anomalous penetration, the 
interaction potential should he weak compared to the particle energ 
and should he localized around the lattice sites. It is shown that 
the more localized and weak is the interaction potential, the more 
pronounced is the channeling effect. For charged particles, two ' 
potential models have heen employed, the screened Coulomh and Bom- 
Mayer. The channeling hehaviour is again determined hy the range 
of the potential. When the width of screened Coulomh potential ^ 
is equal to the width of Bom-Mayer potential, the screened Coulomh 
potential is more favourable for channeling. The temperature de- 
pendence of the process due to lattice vibrations, comes primarily 
through the Behye-Waller factor and tends to reduce the anomalous 
effect. The energy dependence of the expressions supports the | 
experimentally established fact that thickness into the 

crystal at which one half of the initially channeled particles , 
have escaped, is proportional to the particle energy. 



Tbe problem of emission of particles from the source embedde 
in the crystal has been studied quantum mechanically in Chapter I?. 
The real part of the renormalization matrix adds a correction term 
to the -B'ourier transform of the potential, which determines the 
width of the emission pattern. The magnitude of this correction 
term is extremely small at high particle energies but starts increa 
ing as V^p ^"0^ energies. Unlike the previous quantum mechani- 
cal treatment of this problem, the attempt here has been to calcu- 
late the attenuation due to inelastic processes using renormaliza- 
tion matrix and is found to be small compared to phenomenologically 
chosen value. The attenuation depends upon the particle mass 
and energy as (m^/Ep)"'/^, in interesting feature of the pre- 
sent formalism is qualitative indication of difference in widths 
of electron and positron emission patterns even in the two beam 
theory. 

The unergy loss of the particles moving in channels is 
mainly through the excitation of the conduction electrons as colle^ 
tive oscillations. The corresponding contribution has been calcu-, 
lated numerically in Chapter IV using a dielectric formulation 
with different models for dielectric function. The stopping power 
dE/dx is foimd to be very nearly inversely proportional to velocit 
for particle velocities less than about ten times the Eermi velo- ' 
city. Eor higher velocities, the- usual velocity dependence | 

(dE/dx oc 1/Y ) is obtained. Einally the conclusions have been | 

summarized in the last chapter. 



CHAPTER I 


INTRODUCTION 


^ Discovery and Applicati ons of Channeling ; 

The study of channeling and blocking phenomena in crystals 
is a new and developing branch of Solid State Physics. When 
energetic particles incident on a crystal move along directions 
close to the principal crystallographic axes or planes, they 
penetrate anomalously long distances and are said to be channeled. 
The importance and applications of channeling technique in Solid 
State Physics is attracting many workers to look into the details 
of the phenomena and to search for new applications of the tech- 
nique. A substantial amount of experimental work has been done 
during the last five years; mostly at the Chalk River Nuclear 
Laboratories, Canada and the University of Aarhus, Denmark. The 
technique is finding extensive use and has proved to be very 
useful in Nuclear Physics too, as for example in nuclear life- 
time measiirement s. 

Historically, the directional effects were anticipated 

-I 

as early as 1912 by Stark and 7/endt but the detailed theoretical 

and experimeutal investigation of the problem was delayed until 

2 

the last decade. In 1965, Robinson and Oen made a computer 
simulation of the heavy ion motion in ciystals. They investigated 
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the slov/ing down of 1-10 keY heavy ions in a lattice model which 

included realistic repulsi-ve interaction and obtained very large 

penetrations for particles with initial velocities nearly parallel 

to close-packed directions. Experimentally this effect was first 

observed in studies on the penetration of keY ion in crystals by 
5 4 

Piercy et al. and by Lutz and Sizemann . They found by measur- 
ing the distribution of penetration distances, that a significant 
fraction of ions incident along low index directions had anoma- 
lously long range. At about the same time, Nelson and Thompson 
observed anomalously high transmission of 50-75 keV protons and 
helium ions in low index directions of thin gold single crystals. 
The effect was fuorther confirmed by the experiments of Deamaley 
and others who reported that in the transmission experiments, 
using thin single crystals of Si, 1-10 MeY protons and alpha 
particles lose energy more slowly, if injected parallel to a 
close-packed axis or plane. Similar reductions in energy loss 
were observed'^ for bromine and iodine ion beams at energies 
upto 100 MeY, These experiments and the detailed range measure- 

Q 

ments of Komelson et al. on heavy ion beams in single crystals 
fully confirmed the channeling concept and established that in 
an aligned crystal, more than 50% of the beam becomes channeled 
and penetrates much further than the observed range for an irre- 
gular array of target atoms (i.e. amorphous target). 

The first accurate and comprehensive theoretical treat- 
ment of these directional effects was given by Lindhard^ in 1965. 
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This trea,tmen,t is based on the classical orbital picture of the 
actual steering mechanism for the motion of the particle in a 
ciystal. The experimental development was extremely rapid follow- 
ing this classical theory of lindhard. 1 good amount of experi- 
mental 'Work has been done on heawy ion ranges in tungsten single 
1 0 

crystal , wide angle Rutherford scattering in tungsten and sili- 

11 12 
con , and other diamond-type lattices (such as diamond, Si, 

Ge, GaP, GaAs, GaSb) , x-ray production in A1 and Cu by proton 
bombardment and neutron production by nuclear reaction of protons 
in copper”^^ (Cu^^ (p, n) Zn^^). A more detailed account of expe- 
rimental work has been reported a.t the International Conference 
on Electromagnetic Isotope Separators and related ion Accelera- 

1 5ci 

tors and their Application to Physics and at the International 
Conference on Atomic Collisions and Penetration Studies with 
energetic (kcV) ion beams held at Chalk Rivor, Onto,rio . 

The phenomenon of ’channeling’ has proved to be a very 

useful tool for making some very important experimental studies. 

It can be used to find out whether a doped impurity has occupied 

16 17 

a substitutional site or an interstitial site ’ in the crystal. 

18 

further, one can study the lattice disorders and surface 
1 ^ 

effects using this technique. It has been applied to the study 

20 

of nuclear-life times also, in a very interesting way. One 
particularly useful application of the channeling technique has 
been the study of ion implantation in semiconductors'^ ’ ’ in 

which the location of the implanted ions and a determination of 



4 


damage produced is studied by means of tbe channeling behaviour 

of MeV projectiles. Ihe technique has also been used in finding 

out the location of inert gas atoms in some crystals (such as 

KCl, Cal’2 UO2). Another related phenomenon, which on one 

hand, gave rise to the feeling of the existence of the directional 

effects in ciystals^^’ and initiated its discoveiy and now, on 

the other hand, uses channeling as a technique for its own detailed 

27 

investigation and development , is the phenomenon of 'sputtering' 
where the yield of the atoms from the target due to bombardment 
of high energy ions and protons depends upon the crystallographic 
directions. A simple theory of 'sputtering' has been given by 

00 

Lehman . A more detailed information regarding the discovery 
and applications of channeling and blocking phenomena are dis- 

pQ 

cussed in the review articles by Datz et al. ^ , Erginsoy, and 
30 

Davie s-^^. 

The following section contains a brief review of the clas- 
sical theory proposed by Lindhard. The underlying basic assump- 
tions which lead to a consistent approximation method for treat- 
ing the channeling phenomena, have been given and the conditions 
of their validity are derived briefly. The applicability of 
this theory and the cases, vdiere a quantum mechanical theory has 
to be used in view of the fact that the former cannot give any 
thing beyond the gross features of the phenomena, are indicated 
in the last section of the present chapter. 
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• 2 Review of the Classica l Theory : 

The classical theory of directional effects in perfect 

q 

crystals proposed by Lindhard is based on the 'orbital picture' 
of fast moving particles and their steering by the densely 
packed strings or planes of atoms existing in the crystal. The 
theory has been constructed on the basis of four main assump- 
tions which lead to a consistent approximation procedure for the 
treatment of motion of the particle in channels. 

Rirst, the angles of scattering of the particle may be 
assumed to be small. This assumption is valid for fast heavy 
ions moving close to the crystallographic directions and obviously, 
scattering by large angles would imply that the original direc- 
tion is completely lost. The deflections can then be computed 
classically as if the atoms are infinitely heavy i.e. as if the 
particles were moving in an external fixed potential. 

The second assumption is that the strongest correlation 
occurs for a collision with a row of atoms. The other atoms in 
the crystal have a negligible influence on the correlations 
occuring in one particular row of atoms. This assimiption holds 
because a collision demands that the particle comes close to the 
atom and collisions are strongly correlated if the particle moves 
at a small angle with a row of atoms (as a consequence of the 
first assumption). However, if it passes close to one atom in the 
roWjrt must also pass close to the neighhouring atoms in the 
same row. This leads to the concept of a string of atoms, 
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characterized merely by the constant distance of seioe.ration 'd* 
of atoms, placed on a straight line. In first approximation, it 
can be said that the collisions occur with one string at a time, 
the string collisions being independent and uncorrelated. The 
physical importance of the string is emphasized by the fact that 
practically all the physical processes caused by the particle or 
influencing its path, demand that it comes close to the string. 

One exception is the resonance excitation of the atomic electrons 
which may take place far av/ay from the particle, if it has a 
high velocity. The simplicity of the string approximation is 
due to the fact that the lattice structure does not enter into 
the problem, the only lattice parameter coming into picture 
being the distance d between atoms in the string. The strings 
belonging to low index directions have a small value of d, and 
are the most pronounced ones. Correlations weaker than those 
of strings are expected for crystal planes, atomic pairs etc. 

The third assumption is regarding the adequacy of the 
classical orbital picture. Since the individual collisions of 
the particles with charge z^e and velocity v with the crystal- 
atoms of charge Z2e, need quantal corrections either when the 
quantity x = (2z^Z2e'/hv) is not large compared to unity or 

when the impact parameter is large,* the classical approximation 
might seem to be doubtful in several cases. However, it has 
been shown by lindhard (Ref. 9 , Appendix B) that the cla,ssical 
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description of many successi-ve collisions with atoms in a string 
does not become invalid at high velocities. 

lastly, the idealized case of a perfect lattice and a 
perfect string may be used as a first approximation. The ther- 
mal and zero point vibrations of atoms can be considered to be 
comparatively small and their influence on the symmetry proper- 
ties of the lattice is assumed to be negligible. 


With these assumptions lindhard introduced the idea of 
continuum approximation to the potential of a string. Using this 
in conjunction with the assumption of small angle scattering 
leads to an expression for critical angle for channeling. The 
basis of the continuum approximation is to assume that many con- 
secutive atoms in the string simultaneously influence the parti- 
cle trajectory. The average potential responsible for the deflec- 
tion is obtained by averaging over the whole string. Thus if we 
take ions as the incident particles, the average potential at a 
distance r from the string is given by 


U(r) 



—OO 


( 1 . 1 ) 


where V(R) is the ion-atom potential and d is the distance between 
atoms in the string which has been assumed to be along the z-axis. 
Eor R not very much larger than a, the potential T(R) is essen- 



V(R) 


■t o (R/a) 


( 1 . 2 ) 
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where a is the screening length of the particle atom interaction 
and is given hy a = a^ x 0.8853 and ^^(R/a) 

■zp 

is the Fermi function belonging to .one isolated atom . Equation 
(1.1) shows that, the variation of U with 1/r is by one power 
less than the variation of 1 with 1/r. From the equations (1.1) 
and (1.2), we may write 


z.z,e 

U(r) = -L±- g (r/a) 


(1.3) 


with 


C (r/a)= / 


dz 


z^-fr 


(j» z^+rVa) 


(1.4) 


A detailed and more accurate estimate of S (r/a) is given 
elsewhere^^. For the present purpose, we may use some what simp- 
ler estimates. In order to get qualitative insight in the beha- 
viour of ?, i.e. of U, we note that 4^ ^ 1 for R/a << 1 so 
that from (1.4), (r/a) should increase logarithmically for 

small r. Thus, 


C (r/a) % 2 log 


Ca 


for r < Ca 


(1.5) 


where the constant of integration 2 log C is determined by the 
screening. A better estimate valid for all r has been taken to be 
be /Ca^2 


C (r/a) ^ log [ (^)^ + 1 ] 


( 1 . 6 ) 


so that ^ (r/a) (Ca/r) for r > Ca. The constant C has been 
chosen to be ^f3 which gives fairly good overall fit. The expres- 
sions for the density of electrons p (R) and atomic potential 
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V(R) corresponding to formula (2.6) for C (r/a) are given by 

,2 


p (R) = ^ ^2 


(Ca)' ^ 


V(R) 




[ R^+(Ca)^ } 


(1.6a) 


(1.6b) 


The condition for validity of continuum approximation for 
the potential due to the string is obtained by demanding that 
the distance d between successive atoms in a string is small 
compared to At. Vj| where Vn = v cos is the velocity component 
parallel to the string and can be taken approximately equal to 
V for small 4^ > and collision time At = ) 5 1 

being the impact parameter with the string and (l),the cor- 

responding minimum distance of approach. The condition for vali- 
dity of the continuum approximation thus becomes 


At. V cos ^ 


r^,- (1) 

min ^ ^ 




> d 


(1.7) 


let us apply the condition (1.7) in its most restrictive form 
so that we demand its fulfilment for 1 = 0 and determine 


r . (1=0,4;) (hereafter written simply as 3? ^„) by 

iHLX«X«Xi, XXaJLXX 


' min'^ 


12 2 2 
sin ^ ^ f 


( 1 . 8 ) 


so that the condition (1.7) becomes 
C| 

V 


§ exp ( - ^ ) > 1 


(1.9) 


P 

with b = z^z^e /E and E is particle energy. For ^ increasing 
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from zero, the inequality is -violated first hy the rapid decrease 
of the e:Kponent ial, provided (Oa/ipd) can remain large. Thus the 
condition (1.9) requires that 

= VWEp = = 2z^Z2eVd (1.10) 

provided Ca/ il^^d > 1 i.e. < a/d or > E' = (1.11) 

At low energies where (1.10) is not valid, one uses (1.6) 
in (1.7) and ( 1 ,8), getting the condition 

^ < ip 2 = (Cai|;^/df2)’’/^ (1.12) 

and since C/f2 % Y3/]F2 _y 1, the critical angle \J; g applies when 

> a/d or E < E* (1.13) 

Erom the above expressions ( 1 .S)-( 1 , 13) , we note that the poten- 

2 

tial energy barrier E (for E > E’) is independent of energy 

ir ' ir 

2 1/2 

but for Ep < E’, the corresponding barrier Ep'1^2 "paries as Ep ^ 
and decreases as the particle energy decreases. In contrast to 
in (1.10), ^2 (1.12) depends upon the atomic radius and 

on the behaviour of screened atomic potential. As a consequence, 
(1.12) cannot be expected to hold accirrately at very low energies. 

The behaviour of a beam of particles moving through a 
lattice is qualitatively discussed by equations (1,10) and (1.12). 

If the initial angle is less than the continuum picture 

of the string applies and the beam obeying this condition is said 
to be aligned beam (the condition for aligned beam being ’I' < f 

O’ of order 1-2) ?dien tp > C'4> ? one gets the random beam. Thus 
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the yield of a reaction requiring the incident particle to come 
close to the lattice atoms (within about 0.2 1 ) such as Ruther- 
ford scattering, nuclear reactions etc. , will depend upon the 
angle of incidence, it starts decreasing when the angle of 
incidence with respect to a prominent ciyst alio graphic direc- 
tion is less than and it reaches minimum walue for if' = 0, 

corresponding to a completely aligned beam. This minimum yield 

■54 

has been estimated to be 

X = Nd-rc (p^ + a^) 

min X 

where is the mean square amplitude of lattice vibrations 

perpendicular to the string and a is the Thomas Fermi screening 
distance . 


Another Important case of channeling is the case of planar 
channeling which occurs when the particle moves parallel to a 
ciystallographic plane. The correlations here are weaker and 
less effective than for the string case. Using continuum appro- 
ximation, the average potential in this case can be written as 

00 na il y w 

Y(y) = Nd f 2nr dr V(f y^+r^) (1.14) 

o 


where y is the distance from the plane and Wd represents the 

ir 

average number of atoms per unit area of the plane, d^ being 
the distance between the consecutive planes. The ion-atom 
potential V(R) is given by (1.2). Using the standard atomic 
potential (1,6b), one gets 
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1 /2 

Y(y) = 2-KZ^z^e^ Ndp [(y^ + C^a^) - y ] (1.15) 

which corresponds to the expression (1.6) for 5 (r/a) for the 
string case. 


The criteria for the use of continuum potential for planar 
case is much more involved than for a string. The approximate 
condition for this case, as found by lindhard^, is 

^2(i^^2)2 j-^^2^^^1/2 _ ^ j > 

with a 


a 


(1.16) 


' , z -z^e 

E /E , C = y . /Ca, E ' = 

p / p» "^miu ’ p 


27tC^a\d. 


z^z^ X 30 eV. 


■p 


and 


^min minim;im distance of approach to the plane. 


Another consequence of directional effects is observed 
when the particles, emitted from sources (such as radioactive 
nuclei) embedded in a crystal propagate through it. The inten- 
sity of the emission pattern is found to depend strongly on the 
direction of observation with respect to a prominent crystallo- 
graphic direction or plane, Eor positrons, the emission pattern 
has a dip along the directions along which maximum channeling 
occurs (when the emitter is at a lattice site). On the other 
hand, one gets a strong peak along these directions, for the 
case of electron emission. This phenomenon is called blocking 
and has also been treated by the classical theory of lindhard. 
Supposing that the emitting nucleus is in the neighbourhood of 
an atomic position in the perfect string, one takes the probabi- 
lity distribution dP(r) of this nucleus, in space to be of the 
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Gaussian type 


ap = 


2rdr. a/p' 


(1.17) 


wTaere is the mean square displacement due to the thermal 

_r2/p2 _1 

vibrations and a=(l-e ° ) 


2 1 

and r^ is defined as ixr^ = (Ud)" 


is normalization constant 
Then the probability dis- 
tribution of transverse energy of particles is found to be 
(assuming P « r i.e. a 1 ) 


^ p2 2 ^ _1 

n [- - 75 - ® - 1 ) l-exp[- ( 1 . 

^ P 


and the integrated dip in the emitted intensity is given by 

2 


Tt 4'. 




1 


log 


rC^a^ + p^ 


, Y= hl8 


(1.19) 


This expression shows that as the thermal vibrations tend to 

2 

increase the mean square displacement p , the dip becomes less 
pronounced. For the ideal case of no theimal vibrations, one 
gets maximum dip. 


Equation (1.19) can be used to calculate ^ if P is known. 
On the other hcmd if Q is experimentally measured then the cal- 
culated value of p from eqn. (1.19) can be compared with its 
value determined otherwise. Apart from such estimates, it is 
seen that is proportional to 'f'.j = 2z^Z2e /dE. Thus the par- 

ticles emitted with less energy will show a more pronounced dip 
than that exhibited by particles emitted with higher energy. 


18) 
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1 . 3 Meed for Quaxitum Theory “ 

The classical theory reviewed in the last section as pro- 
posed hy Lindhard has been found to he in good agreement with 
experiments on heavy ion and alpha particle channeling. The mag- 
nitude of critical angle, ranges of ions, yields of Rutherford 
scattering and nuclear reactions etc. have been measured and 
results compared with the classical predictions are found to he 
in satisfactory agreement “ However, for light particles, 
such as the electrons and positrons, the classical theory cannot 
he used if one needs anything heyond the gross features, mainly 
because of two reasons: Firstly, the phenomena must he closely 
coupled to the Bragg reflections and other related interference 
effects which are observed in electron micrographs. Actually the 
phenomena observed in the electron microscope is described in 
terms of wave interference; Bragg angles and resonance widths 
dominate the intensity patterns even when the electron wave-length 
is veiy much smaller than the lattice spacing of the target ciys- 
tal. Secondly, according to classical mechanics, the trajectory 
of a particle with a definite energy, moving through a fixed po- 
tential field ■V(r) remains the same if the mass of the paarticle 
is changed. In the corresponding quantum mechanical treatment, 
there is a direct dependence on the mass which vanishes only 
when the particle mass becomes arbitrarily large. It has been 
demonstrated by lervig et al.^^ that electrons and positrons 
definitely exhibit many of the non classical features. 
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Actually, it has heen established by DeY^ames et 
that in the phenomena of channeling, the particle mass and the 
strength of the interaction potential simultaneously play an 
important role in determining the extent to which the classical 
theory may be used. They have shown explicitly that when, either 
the particle mass or the strength of the interaction potential 
becomes small, one should use quantum mechanical treatment. This 
criteria indicates clearly as to why neutrons and protons, ins- 
pite of their equal mass, aie to be treated differently. Because 
of the weak strength of interaction for neutron case, it has to 
be treated quantum mechanically where as the protons, which inter- 
act via strong Coulomb (or screened Coulomb) field may be treated 
classically. Similarly, the electrons and positrons are expected 
to exhibit pronounced quantum effects in the emmission expert- 
ments"^ ’ because of their small mass. Even in the case of 
protons, where the classical treatment should be applied and the 
wave mechanical treatment has been shown^^ to yield classical 

results in the limit of proton channeling conditions; there are 

41 4? 

some experimental observations ’ which can be explained only 

on the basis of wave mechanical theory. This point has been 

45 

emphasized by Chadderton by considering the proton channeling 

41 

experiments of Gibb son et al. . The presence of anomalously 
high energy loss and anomalously lov/ energy loss components in 
the transmitted energy spectra is a definite indication of the 
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importance of interference effects due to wa've nature of protons. 

Another evidence for existence of proton waves, is observation of 
42 

star patterns . Apext from these reasons, a purely quantum me- 
chanical description of the atomic processes leading to directional 
effects in crystals is very much desirable from the standpoint of 
solid state theory. 

The wave interference in electron channeling has been treated 

by many authors^'^““^'^ and the problem has been essentially that of 

electron diffraction. On the other hand, for the case of neutrons, 

the effects of crystal periodicity have been neglected till the 

pioneering work of beWames et al.'^^ in this direction^ published 

in 1966. They used quantum mechanics to explore the possibility 

of anomalous transmission of particles in the crystals. This has 

50 

been further investigated in detail by the author-^ . 

In the present work, a general quantum mechanical formalism 
has been developed in Chapter II which has been applied to the 
channeling of electrons, positrons and neutrons in Chapter III. 

The problem of emission of charged particles (such as electrons 
and positrons) has been treated in Chapter IV and the correspond- 
ing blocking effects investigated. The energy lost to conduc- 
tion electrons by the charged particles during the phenomena of 
channeling and blocking has been calculated using a dielectric 
formulation in Chapter V. The conclusions have been briefly 
summarized in the last Chapter. 



CHAPTER II 


QUANTUM THEORY OE CHAITSELIFG 

2*^ Genera l Eormalation ; 

The basic concept in the theory of directional effects in 
propagation of particles through crystals is to realize the cru- 
cial role of the crystal sjTnmetry and its effects on the parti- 
cle wave function. It has been customary in the theory of ano- 
malous transmission for particles to draw analogies from the 

corresponding effects for x-rays which is known experimentally 

■51 52 

as • This effect is known theoretically to be 

a consequence of the crystal periodicity which leads to the for- 
mation of standing waves in the ciystal when the Bragg condition 
is satisfied. Those waves which have their antinodes at atomic 
sites are attenuated at an enhanced rate, while those having 
nodes at the atomic sites are negligibly attenuated and hence 
anomalously transmitted. The similar physical picture regarding 
the crystal periodicity is also true for the anomalous transmis- 
sion of particles in the two beam theory. Therefore every atterapt| 
to understand the effect theoretically, must necessarily assume 
the existence of the crystal periodicity and any deviations from 
it must be treated in perturbation approximation. This feature | 

Q 

is present in the lindhard's classical theory^ (discussed in 
Sec. 1.2) in which as a first approximation, regular chains of ; 
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close -packed atoms jcalled strings are assumed to exist in the 

45 

crystal and also in the quantum formalism of DeWames et al. 
where the crystal is assumed to he in ground state ha'ving exact 
periodicity of lattice. 

In the present work, v^e assume that the crystal is ini- 
tially in a -very low lying state ln> , so that the approximate 
crystal periodicity still exists. This effectively means that 
the crystal is at very low temperature and as far as the crystal 
symmetry is concerned, tho lattice vibrations do not destroy it 
to any considerable extent. Therefore, as soon as the particles 
which outside the crystal were described in terms of plane waves, 
enter into the crystal ,they are influenced by the periodic inter- 
action potential of the crystal. It is known from the band theory 
of electrons in solids that as a consequence of Bloch theorem, 
the wave function of a particle in a periodic potential should 
be a Bloch wave. Consequently, the incident plane waves must 
go over to Bloch waves inside the crystal whoso wave vectors 
should be determined by using the boundary condition at the crys- 
tal surface . 

Another approximation that we make, is that only a few 
states which are close to the initial state of the system, are 
excited due to the incoming particles. This may not be obvious 
due to the fact that channeling takes place with fairly energetic 
particles. However, if we realize that the particles moving 
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close to the channeling conditions , do not come very close to 
the lattice atoms, as assumed hy Lindhard to illustrate the 
steering mechanism and hence they have only gentle collisions 
with the lattice atoms. Consequently the approximation is 
quite reasonable under the conditions that exist in a channel- 
ing experiment. Under this approximation, the matrix element 
of the potential between the crystal eigen states |n>and 
|m>has approximate symmetry possessed by V when | m > is closed 
to the initial state |n>and can be expanded in lattice Pourier 
series. In particular, we will consider only one phonon inelas- 
tic processes implying thereby jm> = |n + 1 >. 

llov^' we consider a perfect crystal described by the 
Schrodinger eqn. 

H^in> = ln> (2.1) 

where is the second quantized Hamiltonian and is the 
energy of the crystal when it is in the n-th eigen state, assumed 
to be the initial state which is very low lying at some low tem- 
perature. A beam of particles each of mass m^, energy E^ is 
incident upon it. The total Hamiltonian of the system consist- 
ing of the crystal and the particle can be written as 

H = H^ + Hp + Y (2.2) 

2 2 

where H^ is the free particle Hamiltonian -?i Y /2m^ and Y 
represents the interaction between the particle and the crystal. 
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Thus -the total system is described by the Schrodinger eqn. 

Hf = E'i' (2.3) 

The total wave function 'i' may be expanded as 

'i' = 2 <()™(r) lm> (2.4) 

^ m 

where r is the position of the particle, is the actual 
position of ^ -th nucleus i,e. R , £ and u^ being 

the normal position and the displacement of a-th atom, respec- 
tively. Using eqns. (2. 1 )-(2.4) , we get the eqn. for 4>^(r) as 

Ci^ + (E - E„)] = S V t2.5) 

Eqn. (2.5) is similar to the corresponding equation obtained 
in scattering theory, with right hand side including a non- 
local term in n space i.e. both the elastic (m = n) and in- 
elastic (m ^ n) processes being included. ¥e will see that the 
inelastic part is responsible for both, attenuation and anomal- 
ous transmission. 


As discussed earlier, the particle wave function 4’j^(^) 


should be a Bloch wave and can be written as 

ikT,;r. r 


( 2 . 6 ) 


where the periodic function ■£q(£) crystal periodicity such 

that n^(£) = ■*'?_)> 2 . "being a lattice vector and ^ is the 

wave vector of the particle inside the crystal which may be 
determined by using the boundary conditions at the crystal 

surface. The periodic function n^(£) can be expanded as 
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so that 


Ujj(r) = s 


= 2 u, (n) 
^ h ^ 





(2.7) 


where is a reciprocal lattice vector, low, as indicated in 
the beginning of this section, assnining approximate symmetiy 
for jm> close to jn> , we expand 

iKi^.r 

7 (r) = S V.(n, m) e (2.8) 

h 

Using equations (2.7) and (2,8) in (2.5), we get 


2in [^~h "^n ] ^h^^^ ^ ^ 

0 g,in 

where = 2]!]^ (E-E^)/-?i^ = 2m^ E^ 

energy. 


■7h_g(n,m) Ug(m) = 0 

being the particle 


(2.9) 


Eqn. (2.9) is to be solved for the coefficients u„(m). To 

& 

do this we note that it can be written in matrix form ss 


A A 

A U = -- V U 


(2. 10) 


Here U is a column matrix, each element being specified by the 
crystal eigen stste n and the Eourier index h. Thus, a parti- 
cular element U‘ = u^(n) is h-th Eourier coefficient in the 
expansion of , "the periodic part of the Bloch wave '{'^^(r). 

A A 

Similarly the elemients of the matrices A and V are specified, 

A being a diagonal matrix. The particular matrix elements may 
be written as 

■^gm " 'Sn^ P-h ^vP “ ^n J ^h,g ^n,m 


(2.11a) 
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and 


¥ 


hn 


, 1 , 
Vgt"’'”) = ir 


gm ’n-g^ V 

¥’ "being the •voliime of the crystal. 


(2.11b) 


The matrix equation (2.10) is in general, quite oomplic- 
cated and difficult to sol've. The high dimensionality results 
from the product of the mmber of phonon states excited and the 
nimaber of fourier coefficients to be retained in the expansions 
of u^(r) and • A formal solution of equation (2,10) can 

be written as 


K A 

U = u j_ -- c V u 


■ 1 


( 2 . 12 ) 


where Gr = (A - i£)^_^Q is a diagonal matrix and \]^ is a solu- 
tion of equation (2.10) in the limit of vanishingly weak poten- 
tial, namely when the particle has not interacted writh the crys- 
tal and the crystal is in its initial state ln> , so that a par- 
ticular element of this column matrix can be v\/ritten as 


U 


hm 


u^ (n) 6, 


m,n . 


(2.13) 


The dimensionality of the matrix eqn. (2.12) can be reduced to 
a large extent by using the partitioning technique. This tech- 
nique has been successfully used in investigations of defect 
£4. 

problem ’ . The idea is to partition the perturbation matrices 

such that only submatrix is significant and to keep the rest 
of the elements which are negligible, as null submatrices. Thus 
if we take only one phonon processes, the corresponding sub- 
matrix being represented as v, we can write 


A 

1 


A 

V 


0 


0 0 


(2.14) 
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Here, the matrix elements corresponding to higher phonon pro- 
cesses have been assimaed to he negligible compared to one phonon 
elements contained in Vi Partitioning other matrices accordingly, 

A /-V 

A> S 0 tl Ul * 

& = ( . U = ( - ) and U. = ("";■) (2.15) 

0 G U’ 2i 

where g, u and Uj_ have the same dimensionality as v and corres- 
pond to the same space of n and h as v. Using these partitioned 
matrices in (2.12), we get 


A. A 


and 


u = Ui - g V u 


U' = U^ 


How writing vu = _s and v = s^^^ in (2.16a), we get 

u = - g s 

which yields, s = MUj_ with M = (I + vg)”"’ v 

so that a partictilar matrix element can be v/ritten as 


- u^,(m)- ^ E — -- 2 - ^ 


«r' 

gm ’ 1 

m i l . 


(2. 16a) 
(2.16b) 

(2.17a) 
(2. 17b) 


and using eqn. (2.13), we get for m ^ n 

2m. 






,rh'm 
2 


ivriu“ u„(n) 


S (£h,+ V - ^ 


(2.18) 


The Bloch states (p „(r) excited due to propagation of particles, 
are obtained by using (2.18) and multiplying it with e 
and summing over h', as 

2m 


(j) (r) = - 


h h',g 


M^“ u fn) 


- ISn - 1® 


(2.19) 
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Now writing eqn. (2.9) as 




■g 


h ' 5 m^n 


¥. 


h~h 


(n,m)u^,(ni) = 0 


( 2 . 20 ) 


and substituting for u^,(m) from (2.18), we get 
2 


2m, 


o 


[( 5 - 




- ^n] = 0 (2 


where 


°bg(“) 


2m. 


'B' 


m^n h ’ 


..hn ,.h’m 
^h'm “an 

(K^,+kj,?-k2 - le 


(2.22) 


We now examine the conditions under which only a few wawes 
are needed in the periodic expansion of 4^(r). For the ideal case, 

in which the crystal is assumed to be exactly in the ground state, 

49 

the conditions have been examined for 4^. We shall follow a 
similar procedure. If in the expansion of 4^(r) , only first 
term ■^^(ii) is to be significant then the solution forui^(n) must 


be small compared to UQ(n). Thus 


Uj,(n) 


\{n) + 


ho 






- '43 + 


« 1 (2.23' 


where u^(n) has been written in units of UQ(n) . Clearly, except 
at a Bragg condition, one is dividing a Fourier coefficient of 
the interaction potential (since C^g(3a) << ¥^(n)) by an energy 


of the order of 


\ i: K 

o 
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since s-s discussed iDelow. For reciprocal lattice -vectors 

perpendicular to this is about 50 eV for electrons and about 
0.025 eV for protons and neutrons. Away from the perpendicular, 
but still not satisfying a Bragg condition, for >> K^, the 
second term dominates, tending towards Ep(Z^/k^). Thus at 
sufficiently large energy, only the planes parallel to k^ (so 
that = 0) contribute appreciably to (J)^. To justify the 

one wave picture, when such planes exist, it is necessary that 
the Bourier coefficients of the potential be small compared 
to Bor sufficiently high order reflection 

will in any case be small compared to Y^, limiting thereby the 
number of waves which must be taken into acco-unt. 


V/hen the inequality (2.23) is violated for a. single ref- 
lection - e.g. when a Bragg condition (jZ-j^+k^l = k^ ) is ful- 
filled, the first term in the denominator of (2,25) becomes neg- 
ligible and the corresponding coefficient u^(n) becomes signifi- 
cant. This is the well known two wave picture. There will exist 
two solutions for 4)^ having ^ essentially parallel to k^. 

When such conditions are satisfied, we can writu from eqn. (2.21), 
the two equations for u^(n) and u^(n) as 




,2 ,2 






[|r Ooh(“)]"h(“)=o 


.2 .2 


J 


(2 
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Writing ^ = ^ + 11 where n is small and its squares and higher 
powers will he neglected, we get from eqn. (2. 24) , (using the 
Bragg condition |K^+k^j = and the symmetiy properties'^ of 
the renormalization matrix) , 

Bjj(ii)+0jjPn)]upn)+[2n7j,V’^o(n)+C^pn)Ju^(n) = 0 

(2.25) 

[2nYjjVT„(n)+Coo(n)]u^(n)+[j^(n)+0j,o(ii)]uj,(n) = 0 

where = Si^.n, n being the unit normal to the entrance 

A 

surface and is unit vector along k^. Bor nonvanishing 
values of ■^^(n) and u^(n), the determinant of the equations 
(2.25) should vanish. This gives, on simplification, two 
values of p as 

Thus the two waves have wave vectors (k^+D^) and (k^+I3 )• Ike 
decay (attenuation) of these two solutions with penetration 
distance is determined by the imaginary part of p 

Im p^ = (ky2Y^Ep)[j: Im C^^(n)-Im C^^(n)] (2.27) 

We notice that for Im C. „(n) v Im C^„(n), the solution corres- 
ponding to P will propagate almost without attenuation while 
other is attenuated at twice the rate of one wave solution. These 
solutions are, 

i J?" x* 

c C®^P ~) 

iK.h*£ 


(2.28) 
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Here (r) which corresponds to ? has nodes at the atomic 

sites where the scattering out of the heam occurs, so that the 
attenuation is minimized, while the other wave (Ij" (r) has anti- 
notes at these sites and the attenuation is enhanced because of 

scattering. This picture is analogous to that applied in the 

52 

theory of anomalous transmission of x-rays . It should be 

noted that even when several orders of reflection from the same 

plane contribute appreciably to , these two solutions and 

dominate and one still gets anomalous transmission. In many 

interesting cases, this two beam theory is sufficient to give 

the essential features of the phenomena. A detailed account of 

the many beam solution in electron microscopy has been given 
55 

elsewhere^"^. 

From the discussion following eqn. (2.22) upto this point, 
it is clear that calculations of penetration depths and the con- 
ditions for anomalous transmission need a detailed knowledge of 
the renormalization matrix. However, one must get into numerical 
computations regarding inversion, diagonalization and multipli- 
cation of matrices to obtain further results. This expression 

(2.22) may be seen to be similar to that obtained in the corres- 

56 

ponding exact approach to scattering theory using T-matrix-^ 
where the problem ends up with numerical inversion and diago- 
nalization of complicated T-matrix. Consequently, we shall now 
use first Born approximation to illustrate the analytical 


re suit s 



28 


2*2 ?irst Born »' Approximation : 

Let us consider the matrix M = (I + v. If the 

interaction potential is weak so that the elements of the 
matrix -v are small and the particle energy is sufficient to 
make sure that the elements of the diagonal matrix g are small 

A A 

such that the product vg may he taken to he small compared with 

Av 

the unit matrix I, then we can me.ke a binomial expansion, to get 

M = V-WgT + Vg-VgV - ... . (2.29) 

This expression admits of successive approximations. The use of 
the first term v yields the first Born approximation result. 
Similarly, the second term corresponds to second Born approxi- 
mation and so on. 


Thus the first Born approximation in (2.22) yields 
2m, 




0 


■r^hn ,.h*m 
^h’m gn 


mfrx h ' (K^ , +k^) iS 


2m 


0 1 




2 T 


/dr /dr'e 


h'-M" "m 




m;^n 


h 




where in writing the last step, we have used the expressions for 

hn / \ / s 

Vhijjj etc. from eqn. (2.11) and ( r-r' ) has been suitably added 

and subtracted in the exponential. 

A 

If one retains the second term in the expansion of M, one 
gets the next order correction to as 
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( 1 ) 


Ohg = n Sir /dr'/dr" e " 


X s 2 Y 
m'm^n 


(r) ■^nih (r") 


mn 




X ( 


Wl' 


1^6 \/1^6 


(Kv,.+K.) 2:1.2 _ie 


-)(^ 






) 


^h' '-M^ "-"m v£;i^rr±i^/ 

and in this fashion, the successive approximations may he evaluated. 
While the analytical calculation becomes increasingly difficult 
for higher correction terms, the essential features may be illus- 
trated by first Bom approximation. 


As we have seen, the attenuation of the two waves in the 
two beam theory, is determined by Im n ^ which in turn needs the 
detailed calculation of imaginary part of the renormalization 
matrix elements. From eqn. (2,30) we write 


2m 


Im (n) = - ^ /^ir/dr' e 




. r ' 


"hg 




1 


E (^T S e 
m^ h ' 






8 

(2.31) 


Noting that the imaginary part of ^.r is small over many unit 
cells, we put ^ v. in the expression (2.31). further, since 
the summation over h' in eqn. (2.31) extends over whole of the 
reciprocal space, we can replace it by a sum over k = 
to write 

1 „ i(Kh.+^) •(£-£’) ,2 ,. 2 - 


E e 

h' 


->Sb1 


1 

Y' 


ik.(r-r') p o 
S e 6(k:2-kj^) 

k 


(2.32) 
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Since the wave vector may have any value depending upon 
the incident particle energy, the variable k is continuous 
like the usual wave vector in the reciprocal space, and one 


can replace the summation over k by an integration to obtain, 


Im C, 


2m^% 


(n)5:~— -w- — /dk/dr/dr'e 
V'h^( 27 t)^ 




.2 .2 


^ ''n'n (£’) (2-35) 
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which is exactly the same result as has been obtained by 
stairbing with the first Bom approximation solution to ‘J’ j^(£) 


ik^ r-r' 
m — 


4, (r) = 5 / dr' 

“ 4Ttir-r'l ™ ^ 


in eq.n. (2.5) and then using the appropriate expansions of 

' i ^( e ) atia 

Now let us choose a general interaction potential, which 
interms of the atomic potentials, may be written as "VCr) = 

Z '^a(^!“£a^ representing the interaction potential 

between the incident particle and the atom (or ion) at position 


o. The actual position of the a -th atom may be written 
as Ra ~ 2. E a ’ where u^^ is the displacement from the equi- 
librium position of the lattice site a , due to thermal vibra- 


tions. Thus from (2.33)} we get 


Im 


2m 


P -- /<3k E S V (k-k„-E.)e 

V'h2(2Tx)^ n^n 0,5' ^ 


i(k-kn-K^) • a 


® " 6(k2-k2,)<n|e “ ^ > 


X <n' e “ ^ n> 


(2.34) 
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where la (K) = / Va(r) 

iK.U g 

The matrix elements <n ] e |n’> for n'fn, as they are 

involved in the expression (2.34), correspond to exchange of 
phonons i.e. inelastic processes. The procedure to evaluate 
these, is to expand u^ as. 




■K‘ 

2m 


1/2 


£ 


■s 




d 




if.. 


+ at 


-if. . cr 

-3 -J 


(2.35) 


where and f . are energy and wave vector of j-th phonon, res- 
pectively, and £ ( s = 1,2,3) are their polarization vectors. The 

lattice is supposed to he monoat omic, M denoting the mass of 
each atom. There are N atoms in the ciystal. Here a. and at are 

J u 

the well known phonon creation and annihilation operators for 
3 -th phonon. 

The phonon state |n> is written as jn^jng,... >, n^,n 2 ,... 
representing occupation numbers corresponding to different normal 
modes so that the matrix element to be evaluated may be written as, 


iK.u 

i(Q.-a.+QT at) 

, ...je 3 OCT 3'[ 


<n' 1 e n> = <njj , n^ 

• * * 

h 

i(Q. a.+Q* at) 

je 'V 0 ^OCT o'jn.> 

‘ ' 0 


= n 

0=1 ^ 

(2.36) 

^2 

«3o = ( wr 

J 

1/2 if,.a 

) 2 (K.e )e 3 

s ^ 

(2.37) 


for an Inelastic process involving multiphonon exchange, it has 
been shown by Kothari and Singwi"^' that the probability of any 
two phonons to have same wave vector (i.e. exchange of any two 
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phonons corresponding to the same mode of vibration) is negli- 
gibly small compared to the process in which all the phonons 
being exchanged, have different wave vectors. Therefore, for 
an 1 phonon exchange process, eqn. (2.36) may be written as 


iK. u 


<n’ 


a 


n> 


1 i(Q^^a. + 

n <n . + 1}© ^ ^ ^ ^|n.> 

j=1 3 3 


n 

3 = 1+1 


<n . I e 

J 


i(Q. a.+Q* a^) 

3°^ ^ ^ |n.> (2.38) 

D 


Howe-ver, at any finite temperature T, since the exact n. corres- 

j 

ponding to a given initial state |n> is never Imown, we should 
take an average over initial state so that the required matrix 
element becomes 

iK.U oo oo oo i\f 

<n'|e jn> = £ S . . . £ ( n w(n.)) 

ni=0 n2= 0 0 3 = 1 ^ 


1 

n 

d = 1 


<n .+1 

J 


i(Q. 




2o 2 3a 3 


t) 


1 n . > 
' 3 


n 


<n. 


3=1+1 


KQ, 


.a.+Q. 


at) 


3a 3 3CT 3 


in.>(2 


where co(n.) is the probability of n. phonons of wave vector f. 

J 3 ““D 

being present in the field when the lattice is at a temperature 
T and is given by 


" (“j) 


(1-e 






(2.40) 
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I'or an 1 phonon process this averaging has been shown to be 
equivalent to replacing the operator of the form exp (U) by 
(TJ /l! ) exp (^ <U >Qi). Thus for emission of 1 phononswe get 


<n* I exp 


1 

( iK. u ) j n> = n 

- -a j=i 




-W 


1! 


(2.41) 
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where the Debye-Tfaller factor 2W is given hy^ ' , 

2W = 2DK^ 

2 1 (2.42) 

v^ith 2D = (hV2i\®) 2 ~ coth ( 5 ./2k^T) 

We note that eqn. (2.41) gives the matrix element for 

exchange of 1 phonons each having different frequency. But, 

since it has been shown ^ that exchange of two phonons of same 
-1 

type is N times less probable compared to one phonon process, 
eqn. (2.41) gives the matrix element for a general 1-phonon 
process within the approximation of neglecting terms of order 

_-i 

h compared to unity. 


We also see from eqn. (2.41) that an 1-phonon process 

involves the 1-th power of Q. 'v (energy transfer/il 

Do “ D 

which is a small quantity, so that the corresponding matrix 
element decreases rapidly as 1 increases. Moreover, the multi- 
phonon processes do not permit even the approximate symmetiy, 
which is so crucial for anomalous transmission. Therefore as a 
first approximation, we take only the one phonon contributions 
in eqn. (2.41) so that we can write 

<n ’ I exp(-iK.u^ , ) j n> = -iQ*» (n^+1 ) "*/^ exp(-DK^) (2.43) 

where the subscript j = 1 has been dropped from Q.]^, and?-]. 


The corresponding matrix element for absorption of 1- 
phonons is obtained similarly and one gets 


1 [iQi 

<n|exp(iK.u ) n’> = n - ‘ 

~ -u ' 3=^ 




1/2 




1 


-W 


1 ! 


(2.44) 
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so that in one phonon approximation j we have 


<nl exp (iK.Up)|n'> = iQ^n^”^/^ exp (-DK^) (2.45) 

After substituting such matrix elements in eqn. (2,34) one 
should average the expression over initial state. This is because 
the crystal is at some finite temperature T and due to thermal 
vibrations, n^ , corresponding to a given initial state jn>, is 
never exactly known. Therefore we replace n.j and (n^ + 1) by the 
corresponding averages 


and 


<n^> = 

<n^+1>= 




-1 


(e - -1) 

-5/k^T -1 
(1-e ® ) 


( 2 . 46 ) 


where is the Boltzman constant. Thus finally we get from 
eqn. (2.34), 




im CvgCn) = - T Idk E VgCk-k -L) Vg, (k-t-E ) 

Y'-K^(2-n:)^ o. a* n n n g 


X 4 S( * 


2 ^^6(k2_k2^^) if.(a-o*) 6(k2-k2_^) -if.(£-0' 

( 1-0 ® ) (e ^ -1) 




X exp [i(k-k^ ) . (a_ a ' )} (k-k^-K^ ) . (k-k^-K^) 
X exp(_D[(k-k^-K^)2 + (k-k^-Kg)^]) 
since exp (+ iK-j^.c) = 1 . 


(2.47) 


Since any of the phonon modes could be excited due to the 
propagation of particle through the crystal, we have summed 
over all possible values of f and divided by N, the total num- 
ber of possible phonon modes. The probability factors for the 
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phonon to have corresponding energy C are supplied hy (1-e 


C/1CbT, -■ 


and (e 


5/V' 


) 


• 1 ) 


-1 


Thus eqn. (2.47) gives the expression for 


Im C^g(n) in one phonon approximation, for a general interaction 


potential between the particle and the crystei lattice atoms and 
for a general dispersion for the lattice vibrations. 


In order to proceed a little further with eqn. (2.47), we 
assume the crystal to have no isotopes present and that the nuc- 
lear spin of lattice ions may be neglected. Under this assump- 
tion, the Courier transform of the potential, V ^(K) becomes inde- 
pendent of cr so that writing (K) = 'V’(K) in eqn. (2.47), we got 




/awck-vSh)v{k-vSg)(?E-Sn-Sh)-(^-s„- 


X 


-D[(lc-VSh)2+(k-k^-g )2] l(k-kp.(a-a') ^ 


I, 


Q ry 2in ^ 

~T 

X T-TTr'-ff! e 


o 


T 


2 


? p 2m ? 

) if.(a-a') Ul-) 


"(lie 


+ 


) 




4f . (a -o' ) 


]( 2 . 


If there are Isotopes, however, additiona,! terms will appear in : 
eqn. (2.48). lor example, an isotope sitting at the origin, whose 
interaction potential is oV in addition to that of a normal atom, 
will contribute terms involving integrations over the products 
of Y and 6Y in eqn. (2.48). ! 

How the summation overCT_ and a’ in eqn. (2.48) can be writtei 

, * 2 ., 2 

in the form js e | . This type of summations frequently 

£ 

occur in the neutron scattering cross section calculations and 
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can be performed to yield 


'c b 


2 6(k-k^f+K;^.) 


where v is -volme of the unit cell. Using these simplifica^ 

O 


tions, eqn. (2.48) can be written as. 


2m Tt 


-D[(k-k -K. )^+(k-k -K„)^l s3 
X e ^ ^ - n g 1 ?|.)— I [G^(f) 6(k~k^+f+Ej^, ) 


+G2(f) 2^6(k-k^-f+K^„)] 


(2.49) 


where G-^(f) 




(1-e- 


^2 6(k2-k2-2m ?/h2) 

^2^-^ 2MUC ^ 


(2.50) 


Since energy C (f) is periodic in the reciprocal space, ?(f)= 5 ( f+K^^] 
and the f-dependence of G--|(f) and G 2 (^) only through C( f) , we 
h&ye G^(f) = G^(f + Kj^,) for all and G 2 (f) = 
all K]^»f Therefore the terms in the square brackets in eqn. 

(2.49)) which are being summed over f , can be written as 

= Sjs G.^(f+K^,)6(k-k^+f+^,)+ s 
= E Z +&2(f')6(k-k -f’)] 

]q T ^ J 

= U I t^i(f)6(k-k^+f)+G2(f)6(k-4-f)l (2.51) 

where f is restricted to the first Brillouin zone. 
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Now substituting (2,51) in eq.n.(2.49), using 6-functions, 
+ f) for k- integration and replacing summation over f by 
an integration Y’/(2 tc)^ f df, we get 


Im a^g(r,)=- 


m^ df V(f+K. )V(f+K^) 

__2 (f+K^).(f+E ) 

2 Mv^( 27 t)^ ? (1-e“^/^®^) ^ ^ 


X e 




6(f -2f.k^+ 

2 jy rr \2 


Zm^i V(f-Kjj)V(f-Kg) 


■n' Jz 








2^0 C 
:n- 


5(f +2f.lc 


)] (2- 


Ibis general expression admits the use of dispersion for any 
realistic model for lattice vibrations and can bo used with any 
interaction potential by using the corresponding Fourier trans- 
form in eqn. (2.52). 

2.3 Debye M odel Oalculations ; 

It is well known that the Debye model represents a reason- 
able approximation to the actual lattice vibrations in the ciys- 
tals, particularly corresponding to the low frequency part of 
the actual phonon dispersion. In this model, it is assumed that^ 
the lattice vibrates as if it were an elastic continuum but the 
frequencies of vibration cannot exceed a certain maximum value, 
chosen to make the total number of modes equal to the total num- 
ber of classical degrees of freedom. This is specially true for 
the long wave length lattice vibrations which cannot see the 
detailed lattice structure. If the modes are classified by 
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thfcir frequencies, the frequency distribution function g(w) is 
proportional to co in elastic medium. Yfe choose a maximum fre- 
quency ^ Q a-'t which this distribution is cut off so that the 
total number of distinct modes equals 3H. The corresponding 
dispersion is simply expressed as ? ( f ) = 'RKf) = ^cf v.here c 
is the -velocity of sound in the medium. Thus the maximum value 
corresponds to a cut off in the frequency given by f^. The 
Brillouin zone, which limits the range of allowed values of f, 
is replaced by a sphere of same volume with radius This 

sphere is called the Debye sphere. Since the Debye sphere is 

_ 

to contain U points at a density V’/8 tc in f-space, we must 
have 


i.e. 



fo = (fill Vi' <,)■'/’ 


(2.53) 


A more detailed review of the Debye theoiy has been given by 
Blachman^^. 


In spite of the approximate nature of the Debye model, one 
gets quite satisfactory results in many interesting cases, when 
Debye model is used for lattice vibrations. Particularly, 
whenever the results depend upon an integration over f in reci- 
procal space, it is found to give very good results. The most 
striking example is that of neutron scattering cross -sect ions. 
The total cross-sections (elastic and inelastic jboth) are found 
to agree with experiments quite well. An excellent review on 
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•thermalization of neutrons and applicability of Debye model 
has been given by Kothari and Singwi^"^, 

Dor the present purpose we see thc,t our expression (2.52) 
contains integrations over the phonon wa.ve vector f. Therefore 
the use of the Debye model is reasonably good and is relatively 
simple to yield analytica.1 results. This is so in view of the 
fact that particles moving near the principal crystallographic 
axes, do not come much close to the lattice atoms and very high 
frequency modes are not excited even when tho particle energy 
is large. 

Now writing C = 'ftef in eqn. (2.52)‘5 we get 


Im C, (n)= 

2Mhcv^(2Tc) 


5 / -T I —yfief/kify 


X e 




2b of V{f-K )V(f-K„) 






^ticf/kgT 
2m cf 


S 6(f^+2f.k^^- ~|— )] (2.54) 


so that Im C (n)*'-' — — _ 

°° 2li/fficv^(27t) 

, 2 2mr,cf ^ 

^ 6(f -2f.k^+ 


^-■hcf/k^T” 


and Im C, ^(n); 


( 1-e 


2mcv(2%) 


/fdf V^(f) 

^ 2m cf 

6(f^+2f.k - — ^— ) 




5-f f E 


-D[(f+K^)2+f2] 


6(f^-2f.k^+ 


2m^cf V(f-K^)V(f) 


.,.2 




X f. (f-t,)e 


P 2m cf 

6(f2+2f.V -F-)] 


(2.56) 
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Now the ,36 expressions can he used to calculate of 

eqn. (2.27) for any interaction potential hy using the corres- 
ponding Fourier transform in the expressions for Im ^ho ^ n) and 
Im (n) . In the next chapter we give the calculations of 
Im Cj^^(n) and Im Cqq(^) ^or neutrons and charged particles to 
find the attenuations of the two waves 4 ^ and 4 ^ and the con- 
ditions for anomalous transmission. 



CHAPTER III 


CALCUIATIOK USUfC DIEEBEEIT POTEETIAE MODELS 


Pseud opotent i al - Heutrons ^ 

The possibility of anomalous neutron transmission was 

realized shortly after this effect h.:-i. been identified and 

S 1 52 

understood for x-rays ’ . However, no theoretical investiga- 

tion was made at that time and the only experiment directed 

toward seeing anomalous neutron transmission yielded a negative 
65 

result , apparently due to crystal imperfections. There has been 
no further attempt to experimentally demonstrate anomalous neutron 
transmission. The main difficulty in such experiments is that neu- 
trons have already high penetrating power because of extremely 
short range nuclear interaction between neutron and atomic nucleus 
and it becomes extremely difficult to notice directional effects 
in neutron transmission. However, theoretically the problem does 
carry some interest and recently, DeYVames et al.^^ have investi- 
gated the theoretical aspects of anomalous neutron penetration in 
perfect crystals. Their theory is similar to that presented in 
this work but they assumed the crystal to be initially in the 
ground state at T = 0 and used the Einstein model for lattice 
vibrations. In what follows, we shall treat this problem using 
the formalism of Chapter II. 
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For neutrons, the two particle interaction potential is 
very local, the range being of the order of nuclear dimensions. 
For wavelengths larger than the nuclear scattering length e’, 
one can replace the actual interaction hy a potential whose 
width is intermediate between the scattering length and the 
wavelength. We shall choose a form which passes directly over 
to the Fermi pseudo-potential^^ as the width vanishes: 



where a is the scattering length of a bound atom in the crystal 
and is related to the scattering length of the individual atom a’ 
by the relation a = As - 0 , v(r) of eqn.( 3 . 1 ) 

goes to Fermi pseudo-potential^"^ 


Thus the total interaction potential between the neutron 


and the crystal can be written as 


•V(r) 


0 71 




Z a e 
o 




(3.2) 


2 2 

so that V 0 (r-R^) = - 

0 TD 






(3.3) 


Here a. a tias dimension of length and becomes scattering length 
• (of a bound atom at the lattice site a ) in the case of Fermi 


pseudo -potential (p" 


0). From eqn. (3,3) we calculate Va(K) a 


Va(K) = I V Jr) e 


iK.r 


dr 

0 % ' 


% e 


-e:V4p^ 


(3.4) 
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In tbe approximation of zero nuclear spin of the crystal atoms 
and no isotopes present, is independent of a and we get, 


^(E) = - a (3.5) 

o 

Using this in eqn. (2.56) we get 

im 6^^U) - 2m^cv^ 'f K_gJhcf7l^ ® 


2m cf (f-K^).f 

X 5(f -2f.k^+ 

2m„cf 
0 


X e 


-D’[(K^-f) +f ] 2 

6(f^+2f.k^- 


’fl 


)] 


where 


D' = D + 1/4p' 


(5.6) 

(3.7) 


In the Dehye model, the maximum cut-off value of frequency 
f^ is always small compared to a reciprocal lattice vector K^. In 
eqn. (3.7), we note that the integration over f will have this cui 
off f|^ as the upper limit. Consequently, if we neglect f compare 
to the reciprocal lattice vector in the Dehye-Waller factor, 
this will not cause any significant error (because the exponent 
has second power of reciprocal lattice vector). This approxima- 
tion is usually made in the calculations of neutron cross sect ion i 
without introducing any significant error. Using this assumption 
in eqn. (3.6) we get, 


im C^„(n) 


2 

SMm^cv^ “ ~F t “_g--fecf/kgT -2f ,k^+2m^cf^ 


g'ficf/kBT_^ 


6(f^+2f.k^-2m^cf/h)] 


(3.8) 


low to do the angular integration, we choose as z-axis, 


so that Kj^.f = 'K|^f(cos 0 cos a + sin 9 sin a sin ) where the 
phonon wave vector f has been written as in polar coordi- 

nates and a is the angle between and k^ (z-axis). The coor- 
dinate system has been chosen such that is written as (K^,a,0). 
Thus eqn. (3.8) may be written as, 


, ^^3 2 “ ^h fo +1 ^ f^+K, fxcos a . 

Im G^o(n) = — — — • f fdf / dx [ — 


6(x- 


f+2m c/h 


TF 


) + 


f -K^fxcos a 


1 


n 


Te^^^BT-l)’ 


f-2m^c/h 

-a° )] (3- 

n n 


2n 


where we have put cos 0 = x and used f sin4) dcj) = 0. low since,foi 

0 


high particle energies, is always small compared to (2k^+2m^c/h) 


we can use the 6-function in angular integration to get, 

-D'K,. 2 




^ 2_— -la 


2Mm_cv„ k 


C(i+ 


K-j^ cos a 


2k, 


n 


) / ° f 2 coth ( -JSI ) df 
0 


m^c K-, cos a o 

+ -.2 ^2 -r 

^ -^O j 


2hk^ 


(3.10) 


The calculation of Im OQQ(n) is relatively simple. Actually 


Im O^^(n) is related to the total inelastic cross section by, 


Im ^^^(n) ^ 


a (i) 
2m_Y’ T 


(3.11) 


(i) 


where cr J is the total inelastic cross section. We may calcu- 


late it directly from eqn. (2.54) without making approximation 
of neglecting f compared to in Debye-Waller exponentials. The 
angular integration of eqn. (2.55) yields, 
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Im = - 


f 

Blin^cv;ic^ 0 


/ ° f2g~2D'f (^S|)df (5.12) 

0 


This expression is more accurate than the corresponding expression 
(3.10) for Im Cj^qC^) because it does not use the approximation of 
neglecting f compared to K^, Such accuracy is needed in calculat- 
ing the attenuation of the other wa've 4> ^ when 4 is anomalously 
transmitted. 

In the low temperature limit, the f- integration can also 
be done easily, using the replacement coth (■hcf/ 21 c 2 T)^ 1 + 2 e~^^^'^^B^. 
In this limit, the Debye-Waller exponent D becomes, 


D = 

•^1 SkgBjjM 


[ 1+4(T/e3^)2j 


(3.15) 


and from eqn. (3.10) wo get 


Im Cv,^(n) 


Tth^a^ e r %^^h 2 

2Iv!im^cv‘, k L ' 2h k^“ ^o 

o c n n 


1 cos tt ■.! 

to" (1 +-V-- ■^f5 0 (3--' + ) 

n 0^ 

with 4)2 = 02.'*' V4p • Here is the Debye temperature defined 
as hcf = k-D0T^ and we have neglected the terms of order higher 

O ±5 JJ 

than (T/0-[->)^. from eqn. (3.14), we can write 


-D’ 

Im C^^(n) = im C^^(n) e [l+ 


cos a ^m^c K^i <^0® “ 


21ik f^d+eiv©!) ) 


Im ©^^(n) appearing in eqn. (3.15) is that obtained from eqn. (3. 14) 
by putting = 0 so that Im Cj^Q(n) and Im ©^^(n) occuring in eqn. ; 
(3.15) contain the same degree of approximation. 
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Now the wave ^ will experience no attenuation and will 

be anomalously transmitted when Im v. = o i.e. when Im a (n) = 

+ ho ^ 

Irn C^^(n). This leads to a condition for anomalous transmission as 


cos a 


D- 




-Bfg(1+6TVey) 


] = [exp (D^ Kjj 2 ) _ 13 ( 3 . 16 ) 


The condition (3.16) must be exactly satisfied for complete ano- 
malous penetration of and when this is so, the attenuation 

of the other wave <))" is deteimined by 2 Im C^Q(n) and may be cal- 
culated using more accurate expression (3.12). Since l.H.S. of 
eq.n. (3.16) is very small for high energy particles, it will be 
satisfied only for lowest possible values of (Of course = 0 
is ideal case in which it is exactly satisfied but physically, 

^h “ 0 does not correspond toany plane or axis in the crystal). 

As the magnitude of increases, R.H.S. of eqn, (3.16) increases 
exponentially with exponent where as L.H.S. increases linearly 
with Thus only low index (principal) crystallographic planes 

and directions may produce anomalous transmission. Btoreover, 
since + 1/4p , the width of the interaction potential 

must be small compared to the mean s(|uare displacement so that 
1/4p may be neglected compared with L^, to write = L^. Thus 
the condition (3.16) requires that the interaction potential must 
be localized in the vicinity of lattice sites. 


When the condition (3.16) is not exactly satisfied, the 
attenuation of the beam 4 ^ is governed by. 
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Tttl 


5„2 „ 3 
o 


4MiDoCVc 


1 

3 


(1+ 6T^)(1„e)n.r] ( 3 . 17 ) 


where , „ 

-E3_K^ cos a 5111^0 cos a 

g = 0 fll" ^ ^ 

2hk^f^(1+6TVeD ) 

and the condition (3.16) corresponds to 6 = 1 . Is discussed above, 
shorter is the range of the potential ( p*" small) i.e. more loca- 
lized is thu interaction potential, more accurately the condi- 
tion (3.16) is satisfied. In the limiting case when 1/4p may 
be neglected compared to the mean squa,re displacement 2D^, the 
attenuation, at short wavelength is independent of the poten- 
tial width (p~ ) and is detennined only by the mean square dis- 
placement . 


] (3.18) 


Ihe attenuation of the other wave 4" near the anomalous 
transmission condition (3.16) is determined by 2 Im C^Q(n). This 
quantity may be calculated more accurately than Ira (n) from 

eqn. (3.12), In the low temperature limit, replacing coth (ticf/2kg[ 
in (5.12), by and integrating over f, we get. 


im 0„^(n)= 



r die 


^ SkjTD’Z 

-2D.' 


(1+2e-eD/h+- 


®T) t ? 

_ 2D f ^ 

(e ^ ^ ° -1) 


4Dj_^8Dj_'' L 


+ (2+ "h^c^/21Cj"^[r'^Dj_) exp D^) 

x(erf[f (2f^) V2+ ^ _erf ] )I ] (3. 

° 2k^Tr2D'V/^ ^2k_Tr2D-.nV^ ^ J 


2 2 

I 


2_2 


2^2 


2kgT(2])|) 

z ■ p 

where erf (z) =(2/fn) f exp (-t ) dt. 

0 


2kgT(2Dj_) 


19 
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In all the above expressions, the limit of contact Fermi 

pseudo-potential is obtained by replacing by D^. This limit 

corresponds to maximum possibility for anomalous transmission 

condition (3.16) to be satisfied for low index planes or axes 

so that the wave 4 ^ is anomalously transmitted with minimum 

attenuation. On the other hand, the limit of a broad interaction 

potential (p~ large) is obtained by replacing V4' p . But 

in this limit, the condition (3.16) can not be satisfied. This 

—2 

is because the l.H. S, of eqn. (3.16) being less than 10 for 
typical values of and of interest, the R.H.S. cannot be 
any where close to that if 1/4 p is so large compared io 

T 2 

fr . This is more clearly seen from Table I where (1-6) 

1 2 

becomes about 10“ for 1/4p = 100 D^. Thus the condition for 

anomalous penetration requires the interaction potential to be 
short range and the attenuation under the channeling conditions 
is determined by the mean square displacement. The expression 
for € and the condition (3.16) for anomalous transmission indi- 
cate that maximum penetration occurs for low index crystallogra- 
phic planes and directions. 

5*2 S creened Coulomb Potentisil-Charged Particles : 

When the charged particles (electrons, positrons, pro- 
tons etc.) propagate through the crystal, three types of transi- 
tions in the crystal become important: excitation of nonlocalized 
electrons, excitation of electron cloud surrounding the lattice 



49 


TalDle I; Variation of (1-6) with width of the Potential 

and temperatiire when corresponds to ^100l 
planes in Gu(fcc) 


P 

(in -units of ^ 



T = 5°K 

T = 15°Z 

T = 25°K 

T = 35°Z 

0.0 

0.8081 

0.8145 

0.8271 

0.8461 

0. 1 

0.8101 

0.8165 

0.8292 

0.8482 

1.0 

1.0091 

1.0170 

1.0328 

1.0564 

2.0 

1.6097 

1.6223 

1.6474 

1.6850 

3.0 

2.6026 

2.6228 

2.6631 

2.7236 

4.0 

3.9759 

4.0065 

4.0676 

4.1593 

5.0 

5.7131 

5.7566 

5.8437 

5.9741 

6.0 

7.7937 

7.8525 

7.9698 

8.1456 

7.0 

10. 1935 

10.2693 

10.4207 

10.6474 

8.0 

12.8850 

12.9793 

13.1677 

13.4494 

9.0 

15.8380 

15.9519 

16.1795 

16. 5191 

10.0 

19.0201 

19.1541 

19.4214 

19.8212 
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ion and the vihrational transitions of the lattice ion. The first 
effect has nothing to do with the particular direction of motion 
of pairbicle and contributes a large, directionally independent 
background to the total attenuation which is treated in Chapter V 
of the present work. The latter two effects, because of their 
loca.li2ation, will shov\r anomalous behaviour at Bragg angles. In 
the present work, we illustrate the effect of localization by 
considering purely vibrational transitions, choosing screened 
Coulomb potential as the interaction between particle and crystal, 
in this section and Born-Mayer potential in the next section. 

The ions of the crystal are sitting at various lattice sites 
Their interaction with any external charged particle is expected 
to be a Coulomb potential. But because of the collective oscilla- 
tions of the conduction electron sea, this potential is screened. 
Thus the interaction between the lattice ions and the incident 
charged particle is no more a long range -Coulomb interaction 
but is a screened Coulomb potential, the strength of the screen- 
ing being determined by the density of conduction ele ctrons.Thus 
the interaction potential between the incident charged particle 
and the crystal ions may be written as, 

o ^-A|r-R_| 

V(r) = z.z.e'^ E ^ -2- (5.20) 

^ l£-Sal 

so that YJx-Ra) = z.z^ e^ — (3.21) 

^ ^ 2 , j 

I. 

and V(j(K) = V(K) = Anz^z^e^/(A^+K^) (3.22) 
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Using ( 3 . 22 ) in eqn. (2.56)5 we get 


Im C^^(n) = - 


2 m^(z^Z 2 


Mtocv. 






, 2 

-2f.k^+ 

■35^^™: 


2in„cf 
o \ 


X 6(f +2f.k 


2„ of -Bt(f-Eh)^+f^] ] 


The integration involved here is difficiilt compared to that in- 
volved in eqn, ( 3 . 7 ) for the neutron case. Again we neglect f 
compared to the reciprocal lattice vector in the Dehyc-Waller 
exponent and also in the denominators so that the integrations 
are simplified to some extent. Again choosing k^ as z-axis and 
system axes rotated such that f = (f, 9 , tj)) and = (K-^, a, 0). 
Then the angular integration may be done in a similar way as in 


the neutron case, to get finally. 


Im Ct^ -(n) = -( 


2Tim^(z^Z2e^)^ e 


I'fficVc k^ +\ ) 


m cK^ cos a f^ 

) C — 2?nr — IT ^ 

n A 


+ ( 1 + 


cos a f^ coth('tIcf/2k]3T) 


(A^+f^) 


df] ( 3 . 24 ) 


In the low temperature limit, coth 


(■hcf/2kBT) 2 1+26““^°^/^®^ and 


f- integration may be done, to get ^ 

2Tcm (z.Zpe^)^e“^A 

Im (n) = -( 2 ^ 

Mhev^ k^(A^+K^^) 


m e K, cos a f^ 

)[ in(1+ -^) 


Kv, cos a on ^ 

(A, 1)1 


( 3 . 2 , 5 ) 




+ (1 + 
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1 ^0 

With T) = V2f„(T/eB)(1-e'®0^^) 


(^)-si(^)=oh(-i^)] (3.26) 

OO r 

where, si(x) = -I ^ at ehd ol{x) = -/ 22|_. at are the well 

known tahulated funotlons®"^. In obtaining the last two terms of 
eqn. (3.26), we Taave extended the upper limit of f- integration 
to “ because the contribution from f greater than f^ is small, 
for the integrals involved there, at low temperatures. 


Now from eqn. (3.25) we can write, 

-Vh" 


Iin C^^Q(n) = Im C^^{n) 


(1+k^Va ) 


D 


cos a 


2 h. 


'n 


+(ib^0Kjj oos a/a hjj p/°(A,T)) lnC1 + ^3 C.27) 

so that the condition for complete anomalous transmission is, 

2 


cos a 


2k 


'n 


[it 


^) e “ A '- ll ( 3 . 28 ) 

tlI'^^^(A ,5) A A 


The attenuation for the wate , which penetrates anomalously 


at the condition (3.28), is governed by, 


^n' “■ 


(i.eShs.rl 

exp [- ft 2 ^ p ^ 


Mhev^ Ep A 


(3.29) 

.2 


where gQ 


e-hS 


2/. 2 


p = 1 +Ki,7A 


[lt ln(U ^)] (3.30) 


'n 


2h k^EQ^(A,l) 


A‘ 


The effect of the screening parameter A on the attenuation 

T •vin-onn <''^25'l Thus for A ^ large compared to 
can be examined using eqn. v3.<^lp;. 

j)-\ so that A^» A^ » 
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cos a 


lui Ojj„(n)-Im C^^(n) = -Im 0^^{n)[l-e (1+ -S-gj )1 (3.31) 

n 


and 


SC _ (1 + cos a/2k^) exp (-D.^ 


P 


n- 


a 


(3.32) 


which is independent of ^ . Thus when the screening is "very strong, 
implying thereby that the potential is highly localized, the condi- 
tion for anomalous transmission is independent of the screening 
parameter. On the other hand if A small compared to arid 


can be taken to be of the order of f^, then, 


Im Cj^^Cn) - Im CQQ(n) = -Im 


1 - 


© 


■h 




K, COS a 


n 


m^c IL cos a In 2 0 

][ (3. 


2Hk f (1- ^ tan 
n 0 ^ f 

o 


~T 

-1 ^ 
A 


33) 


so that 0 


SC 


-“A 


l+Kj^^/A 


K, cos a m e K^, cos a In 2 

[dt )t ; (3.34) 


2hk^f (1- Ji tan-'’ ~) 


showing that weak screening cannot lead to an anomalous effect 
because for small values of A , e.l of eqn.(3.34) cannot be close 

ir 

to unity and hence the condition for anomalous transmission cannot 
be satisfied. 


3.3 Bom-Mayer Potential - Charged Particles: 

Another model potential for the interaction between the 

incident charged particle and the crystal ions which has been fre- 
* 7 1 

quently used as for example in calculation of energy loss of 

charged particles during their propagation through channels, is 
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the Bom-Mayer potential, 

V(r) 

■ 1 


AS e 
a 


-h 


(5.35) 

where p" ' represents the range of interaction. Here, Both the 
constants A and depend upon the charge numbers of the inci- 
dent particle and lattice atoms. For this interaction we have, 

a I 






2s 2 


so that \{l.) = !/■(£) = 8tcA + Y?) 

G 

Substituting this in eqn, ( 2 , 56 ) ^ we get 5 


(5.36) 


8 A^p^ df 

Im C, ^(m) = r - 

ho^ ' Mhcv f 




'2^2‘,„2‘:J2T2‘ 


(^+K,).f 2 m cf (f-K, ).f 

" -4-)- ~ 


X 




.2s 2 


6 (fS 2 f.k - 


(e^/%^_-l ) 

2 m_ cf 


■n 


0 

T 


) 


] 


( 3 . 57 ) 


Igal. we neglect f compared to ruclprcnal lattice weotor Kj, in the 
Debye— Waller exponentials and in the denominators, as in the case 
of screened Coulomb potential. Performing the angular integra- 
tion by choosing Ic as z-axis, we get, 

^ - >2 ^ .-DKh2 


If Oj,„(n) 


StcA^P^ m^e'"""'''-i3 m^c cos a f^ 

^ 2* I'n Pc 


+ (1 + 


cos a 




h 


fQ 5 coth ('hcf/2k-R[E) 
) / ^f^df ^ 

0 


7^+f^)^ 


(3.38) 


In the low temperature limit, coth (■hcf/ 2 kgT) 2 ,l+e“^^^/^®^ 


so that 
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f- integration may be performed to yield 

2 2 2 
S-reA 6^ m e m^c K, cos a f 

1“ C ° *> 


MRov ^ p^k"^) 2 a p2( p2 ^ f 


‘0“ ‘"h 

:277F 


+ ( 1 + 


cos a 

-^-- 


) C‘ ^Po.’^d 


(3.39) 


with (p T) = ± 


^ ^ O . r^ 








R -Rc^ 


p Hc^ 


■Rc 


+ ci(~^) sin(~~) - si(~- 7 ^) cos(p™-) 
B 


kB^' 




'V 


+ 


p^dic 

k“f 


6 Re 6 Rc B Rc 6 Rc T 

[ol(-^)cos(-^)+si(-g^)sta{-j-j)]j (3.40) 


Thus, Im CQ^(n) = 


eAo^ jBM 


mrct^ JSi P 


(3.41) 


and from eqn.(3.39) we can write 

2 

Iin Cj^Q(n) = Im C^^(n) 




00 


{1+K^2/p^2)2 


[1- 


K^cosa 




a kjj (p^,T) p^Cpy+f^) 


] 


(3.42) 


The attenuation of the penetrating component (j) ^ is gi'ven by, 


1 exp[. 




where 6 


BM 

'P 


MR cw Y E a 
c n pc 

2 

^h K cos a 

® 2-5 


(3.43) 


(i+k^7Po )' 


raQC cos a 


2kj^ 

2 


2* k F„“‘(P<,,T)py(Pc^+f„b 


] 


(3.44) 



56 


Thu condition for channeling which requires Im C^^Cn) = 
or = 1 5 is gi-ven hy 


cos a 

’SF"'""” 

n 


+ 


(P 


m^c f ' 
o 0 


c 0 




,1) 


'] = L( 


-h^ 

7i 


im O^o(ii) 


2 

e 1 ^ -i](3. 


Now we note from eqn. (3.4-4), that for higher particle ener- 
gies, the second and third term may he neglected and variation of 
the channeling parameter (“l-Sp ^'*) with the range of the interaction 
is essentially given hy, 

^ A (3.46) 

The corresponding expression (using similar approximation) for 
the cp.se of sci^ecned Ooulomh potential is obtained from eqn. (3.30) 
and is given hy, 

-Jj E ^ 

(1-€p^^) = 1 - e ^ ^ /(1+ a2) (3.47) 

The variation of (1-Sp®'"^b with and that of (1-6^^^) with A""' 

at low temperatures has heen given in Tables II and III respectively 


3 . 4 Discussion : 

SC 

If wo closely examine the expressions for (l-O- ) and 

P 

(l-£ ®^'*') and the corresponding channeling conditions (3.28) and 

^ 11 
(3.45), we notice that the range of the interaction (p^" and A~ ) 

must he small for channeling conditions to he satisfied. Appearance 

p 

of in exponential as well as multiplicative factor, at R.H. S. 
of eqn. (3.28) and (3.45) requires that the magnitude of should 
ho small i.e. the particles should he incident along principal 


45 



57 


Table II: Variation of (1-6^^) with width of the Screened 
Conlomb potential A “ and temperature when 
corresponds to {100} planes in Cu (fee). 


A-' 

in units of 

(2D^) 

(1 - 

eSCx . 

units of 10” 

■ 2 

T = 5°K 

T = i5°i<: 

T = 25°E 

T = 35°K 

0.1 

0.8242 

0.8307 

0.8436 

0.8629 

CM 

• 

o 

0.8725 

0.8793 

0.8929 

0.9134 

0..3 

0.9528 

0.9602 

0.9751 

0.9975 

0.4 

1.065 

1.073 

1.090 

1.115 

0.5 

1.209 

1.218 

1.237 

1.266 

0.6 

1.384 

1.395 

1.417 

1.449 

0.7 

1.591 

1.603 

1.628 

1.665 

0.8 

1.828 

1.842 

1.870 

1.913 

0.9 

2.095 

2.111 

2.144 

2.192 

1.0 

2.392 

2.411 

2.448 

2.503 

2.0 

6.855 

6.905 

7.006 

7.157 

3.0 

13.449 

13.541 

13.725 

13.999 

4.0 

21.255 

21.387 

21.650 

22.042 

5.0 

29.436 

29.601 

29.927 

30.410 

6 . 0 

37.388 

37.572 

37.939 

38.481 

7.0 

44.746 

44.941 

45.327 

45.896 

8.0 

51.344 

51.541 

51.930 

52.503 

9.0 

57.143 

57.336 

57.717 

58.276 

10.0 

62. 182 

62.367 

62.632 

63.267 
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Tatle III; 



(in units of 


0. 1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.8 
0.9 
1.0 
2.0 

3.0 

4.0 

5.0 

6.0 

7.0 

8.0 

9.0 

10.0 


Variation of (1-6™) ?/ith tbe width of Bom 
Mayer Potential (p^~"*), corresponding to ^OO^ 
planes in Cu(fcc) 


(1 _ e®') in units of 10 ^ 

T = 5°K T = 15°K T = 25^K T = 35°K 


0.8403 

0.9368 

1.097 

1.321 

1.608 

1.957 

2.367 

2.837 

3.365 

3.951 

12.532 

24.480 

37.487 

49.802 

60.478 

69.221 

76.133 

81.483 

85.581 


0.8469 

0.9441 

1.106 

1.332 

1.621 

1.972 

2.385 

2.869 

3.391 

3.981 

12.622 

24.655 

37.693 

50.032 

60.708 

69.436 

76.324 

81.649 

85.721 


0.8600 
0.9588 
1.123 
1.352 
1 . 646 
2.003 
2.422 
2.903 
3.443 
4.042 
12.799 
24.945 
38. 101 
50.488 
61.163 
69.859 
76.699 
81.973 
85.995 


0.8797 
0.9807 
1.149 
1 /383 
1.683 
2.048 
2.477 
2.968 
3.520 
4.132 
13.066 
25.407 
38.707 
51.159 
61.831 
70.478 
77.247 
82.443 
86.392 
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crystallographic planes or azes. These conclusions are siinilar 
to those derived in Section 3.1, from consideration of anonalous 
neutron transmission using lemi pseudo-potential. A comparison 
hetweon the results for the two potential models used for charged 
particles, shows that the screened Coulomh potential is more fa- 
vourable for chiinneling than the Born-Mayer potential. This is 
seen clearly from Tables II and III where (l-^-p ) is smaller 

“DTWr ^ ^ 

than (1-G„ ) for equal values of A “ and B ~ , showing that the 

penetration depth in the former case is greater than in later 

case. Another significant difference betv/een the two models is 

■1 

that the range A “ of screened Coulomb potential depends'only on the 

density of conduction electrons whereas the parameters of Bom- 

Mayer interaction, A and B , depend upon the charge numbers of 

both, the incident particle and the target crystal. In particular, 

72 

in the Brinkman's version of the Bom-Mayer potential, one has, 

A = 2.58 X 10“^ ^eff^ 

/ 1.5 ao (3.48) 

where = ( 2 ^ 22 )*'^^^, a^ is the Bohr radius and 2^,22 are the 

charge numbers of the incident particle and the crystal atom 

_-i 

respectively. Thus as or Z 2 increases, decreases, result- 

ing in larger penetration depth. Similarly for small values of 
2 ^ and 22 , one gets more attenuation and smaller penetration depth. 
But in the screened Coulomb potential, there is no such dependence 
on the charge number 2 ^ and 22 and the penetration depth is governed 
only by the density of conduction electrons through the screening 
properties. 
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There are few interesting features of the channeling 
phenomena which may he seen to emerge as a consequence of the 
preceding quantum mechanical treatment of the problem. Pirst of 
all, we note that the existence of anomalous effect has nothing 
directly to do with the wavelength of the incident particle but 
is rather determined by two criteria: The interaction potential 
must be weak compared to the particle energy as seen from the 
condition (2.23) for validity of two wave picture and that the 
interaction must be localized in the vicinity of lattice sites. 

The conditions for channeling corresponding to all the three 
potential models require that the magnitude of the reciprocal 
lattice vector K^, corresponding to the lattice plane or axis 
with which the Bragg diffraction is taking place, should be 
small. Thus appreciable channeling will occur only when the inci- 
dent beam of particles moves in a direction close to the principal 
crystallographic directions or planes. The energy and temperature 
dependence of the expressions for the channeling parameter (1-6) 
is only through small factors and as such the dependence shows 
that high particle energy and low temperature is favourable for 
channeling. The energy dependence of the factors multi- 
plying (1-6) in the exponent of ^ agreement with the 

experimentally concluded"^^ fact that the thickness into the 

crystal at which one half of the initially channeled particles 
have escaped from the channel, is proportional to the particle 


energy. 



CHAPTER IV 


EMISSION OE CHARGED PARTICLES EROM CRYSTALS 

4- • Int roduction: 

It was pointed out in Sec. 3.1 that Lindhard’s classical 
treatment is applicable only to heavy ion channeling. Eor light 
particles (such as electrons, positrons etc.) the interference 
effects are significant, irrespective of their wavelength and 
then, one has to use the quantum theoiy developed in last two 
chapters. The experimental situation for light -particle channel- 
ing seems to be poor. Some attempts towards such studies have 

74- 

been made by Wolf et.al., who studied electron channeling in Si 

single crystals by scanning electron microscopy. Similar techni- 

que has been used by Joy et al. ' to study the electron channeling 

in ferromagnetic crystals. Some experimental information on high 

energy electron transmission has been reported by Steeds and Valdre 

77 

More recently, Andersen et al. have investigated the positron and 

electron channeling by means of Rutherford scattering. Inspite of 

the poor experimental situation, theoretical calculations'^® 

55 

using standard many beam dynamical theory , have been made. A 

HiO 

recent communication ^ also contains such calculations for elec- 
tron channeling at high energies in gold, 

■ The lack of availability of experimental data on light 
charged particle channeling wo\ild have forced one to predict 
theoretical results without comparison. However, the charged 
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particle emission problems which are governed by the same phy- 
sical mechanism as the channeling problems, have been interest— 
ing experimentally as well as theoretically. A detailed expe- 
rimental investigation of a-particle emission from the radioactive 
nuclei embedded in the crystals, has been carried out by Domeij®”^’®^ 

O’T 

and coworkers who found dips in emitted intensity in low index 
directions. Gemmell and Holland®^ have reported some measurements 
on the charged particle (proton, Deuteron and alpha particles) emis- 
sion from the crystals. These experiments have been essentially 
explained by the predictions of lindhard's classical treatment. A 
significant work on emission of electrons and positrons from 

the crystals has been reported. by Uggerhoj-'^ and coworkers Re- 

86 

contly Walker et al. have studied the effects of Bremsst rah lung 
radiations on the channeling of relativistic electrons and posi- 
trons. Although, Lindhard's classical theory has been used to 

'20 "ZQ 

interpret the patterns observed by Uggerhoj and Andersen'^ 

in view of the discussion given in Sec. 1.3, it is desirable to 

develop a quantum-mechanical theory taking Bragg diffractions and 

36 37 

the interference effects into account. LeWames et al, ’ have 
given a wave mechanical theopy for this problem, but their treat- 
ment, in addition to assuming the crystal to be an Einstein crys- 
tal at zero temperature, starts with an approximate form (one 
particle approximation) of the many particle Schrodinger equation 
and uses the periodic potential field of the classical model, This 
one particle Schrodinger equation conserves the particle energy 
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"the propagation of the emitted particle through the crystal and 
its subsequent escape from the crystal. We shall also neglect 
the interactions (if any) of the emitted particle with the emitter 
itself and thus we are ignoring the processes that might occur in 
the immediate ■vicinity of the emitter. Under these assumptions, 
we are faced with the problem of propagation of particles through 
the periodic potential field of the crystal. The ciystal is again 
assumed to be in a low lying state jn> at some low temperature, 
as in Chapter II. Therefore our system is again described by 
the Schrodinger equation (2.3). The total wave function may 
again be expanded as in eqn. (2,4)}and the equation satisfied by 
the Bloch states representing the particle wave function, is 


j^2y 2 

r _ ___ — +. ’v (r) — E 1 <}) (r) = - E V (r) <1; (r) (4.1) 

L 2m^ nn^ ^ p J ^n^ ^ m''-^ ^ ^ 


where E^ = E - is the energy with which particle has been emitted 
from the crystal. 


Now to find out the intensity at a point r outside the crys- 
tal we make use of the reciprocity relation^ ’ which states 

that the intensity at a point r outside the crystal due to emission 
at a point r inside the crystal, is equal to the intensity at the 
emitter site when the emitter has been placed at the observation 
point r. Thus 

‘f’n £e^ "" ^-e’ 

which holds as long as one neglects the reflections at the crystal 
surface. If the R.H.S. of eqn. (4.1), representing the 
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renormalization is neglected (because its magnitude is sm.all, as 
we will sec later) the relation (4.2) follows immediately from 
eqn. (4.1). In the general case, this relation may be derived 
by using time reversal invariance on the total wave function 
and using the eqns. (2.1) to (2.5). The reversibility rule has 
been experimentally demonstrated^^ to hold in the emission 
and transmission problems. Therefore, the problem has now been 
reduced to that of finding the intensity at the emitter site when 
the particles are coming from a far away point r. The classical 
treatment given by Lindhard, also uses this reciprocity relation 


to relate the two problems of emission and transmission. The prob- 
lem of particle penetration has been treated in deta.il in Chapters 
II and III. Here, we start; by using the first Bom approximation 


result for I namely, 


<P (r) 




ik T— r’ 


471: r^r ’ 




on R.H.S. of eqn. (4.1) and use the usual expansions (2.6) to (2.8) 
for 4j^(r) and V^^(r) to get. 



V = Y „(n,n). Eqn. (4.4) may be written simply, as 
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(26+ Si,)u^(n) + 2 Ug(n) = 0 (4.6) 

whore 26E = k?, - B 

P ^ o M p 

(*"/a»o) [K/+2£i,-y = SfeEj, (4.7) 

\-g = 

how we use the boundary conditions at the surface, to calculate 

restrict the sums in eqn.(4.6) to IT terms, there will 
be IT independent solutions for 6 and hence for the Bloch- coeffi- 
cients u^; all belonging to the same energy B^. These solutions 
must be superposed to obtain a wave function which matches pro- 
perly on to the plane wave which is incident on the crystal surface 
from the direction of the point of observation. Neglecting the 


reflected wave, one can write the boundary condition as, 

i4.4r lK,..r ^ ySn-i 


^ (r) = E A. e 
3 ^ 


/ . N, iK, . r 
S e - 

h 


(4.8) 


for all r such that n.r = 0 (i.c. at the surface). One satisfies 

qp 

these boundary conditions by choosing , 


and taking A^ = where the bar over denotes complex 

conjugation. In fact, the second term on R.H.S. in eqn.(4.9) is 


just n ^ introduced in Chapter II, for the two beam case i.e. 
when single Bragg reflection is considered so that j has only 


two values. 
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The intensity at the emitter site can thus he written as, 
2 




_ (i) . iiCj^.r 2 

S u ^ ^ e S u-^ ^ e 

D h 

£ ju 2 ^h^^^ ® ”j^+2Re S 

D h " r >D“ 


[up)zu(3)e^"^--3 e-- 

n 




n 


. .(4. 10) 


In the last step, the intensity contributed by the cross terms 
between different eigen solutions has been written separately in 
order to display their dependence on the distance n.r, of the 
emitter from the crystal surface. Their evaluation, while straight 
forward, is tedious and adds nothing interesting to the basic cha- 
racter of the solution. Consequently, in the usual many beam cal- 
37 

culations , a simple geometry for the crystal and for the incident 
beam is assumed which assures that these cross terms oscillate 
rapidly with emitter depth and a collection of emitters spread 
over a small thickness will yield an average intensity which is 
just the sum of diagonal terras in eqn. (4.10). This point becomes 
more clear from the two beam theory presented below. 


The intensity pattern which one observes out side a crystal, 
duo to a source placed inside, exhibits considerable structure and 
one should retain maximum number of terms in eqn. (4.10). However, 
almost all the features which appear, can be understood in terms 
of the behaviour of the intensity which one gets for the case of 
single Bragg reflection (two beam theory). 
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Specializing to the case of an emitter located at a lattice 
site (when emitter is at an interstitial site, the intensity maxi- 
muni becomes minimum and "vice versa), the two wave solution takes 
the form, 


£ 

'J>n('fc)|^ = exp {-[1 ] t/^^" } 

■\.n+y2 ° 

6 

+ (1-x)^ exp { -[1+ — ■ 1 t/ p ” } 
+ 2x(1-x) e ° cos ['',ri+y^ t/ 


where 


-Y [1 

. -in-i 

1— w J. 






fi-^y 



y = 


(Re = 

Ohol/Bp> 

= 

h 


/ Ylf »' 

' 0 

;(Im = 

II 


n = 

(is, 

n ^ 

’ 5h" 


Wjj")-'' 


(4.11) 


(4.12) 


and t is the distance of the emitting atom from the crystal 
surface. In deriving eqh. (4.11), we have made the assumption 
that = Im is small compared to vi^hich is justi- 

fied since Im is quite small compared to the Pourier coeffi- 
cient of the potential 1-^, as we shall see in Sec. (4.4). 

Thus the calculation of j4^(t)i^, the intensity at the 
emitter site (and hence at the observation point) is essentially 
that of renormalization matrix and the Pourier coefficient of 
the potential V^. Im has been calculated in Chapter III. 
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Vie calculate now, Re using same approximations, namely 
one phonon approximation for inelastic processes and Dehye 
model for lattice vibrations. 


^*5 Renormalization Matrix Elements: 


From eqn. (4.5) we can v/rite for the real part of the 
renormalization matrix as, 

2m 

^ T " — ^ ^ dk/dr/dr’ exp[i(k-k -K^) .r 

¥'h^( 27 i;)^ - - - L „ _n -h -- 

(4-.i5) 


k -1^, 


where P(1/x) represents the principal value. In writing (4.13), 
we have used the momentum representation of the free space 
Green function, 


e 

I r-r ' 


(2I)? 


/ 


(k^-k^,-iS) 


lor a general interaction potential V(r)=z (r-E^j), this 

a 

expression may be calculated using one phonon approximation for 
inelastic processes and assuming the crystal to have no isotopes 
and zero nuclear spin, as in the calculation of Im C^g(n) in 
Chapter II, we get, 


Re 


^0 

Mv^(27i;T^ 


-D[ff+K^)^+(f+Kg)^] 


df fV(f+K, )Y(f+K^) 


■S(^ ) + 

f^-2f.k^+2m^C/b 


V(f-Kh)V(f-Kg) 




-D[(f-Kj^)^+(f_K )2] 


R(- 


f‘^+2f.k^_.2m^ 5/^ 


2)K4.14) 
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where we have summed o^ver all the phonon wawe vectors f, as in 
the case of Im C^g(n) . 

Now using the Debye model for lattice vibrations, we put 
C = hcf in eqn. (4.14-). Then, we use the screened Coulomb 
potential as the interaction between the emitted charged particle 
and the lattice ions as in Sec. 5.2, to get, 




TcMhcv^ 




(1-e- 


i— ' 

) f^-2f.k^+2mQcf/ti' 


[a f+K ^) [a ^+( f+Kg ) ^ ) 

^ exp(-D[(f-K^)^^(f-Kg)2]) 

■ '?+2f.k^-2m^Gf/h [a 2+( f_K [A ^+( M J 


(4.15) 


Since the maximum value of f in Debye model, small compared 

to the reciprocal lattice vector K^, we neglect f in the Debye - 
?/aller exponent and in the denominators, to get, 

2m p’(z.Zpe^)^ ^ df 


1 1_ , 

?l2f.k^+2m^cf/h (A^+f^ 


?+2f.k^ - 2mQCf/h (A^+f^) 


(4.16) 


where p ' = 1/v is the number density of atoms in the crystal. 

c 

Evaluating the angular integration, choosing as z^axis as in 
Chapter III, we get, 
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Re Ojj^(n) 


Itto (a'®+Kjj2) 0 (A 


^icf 


K^cos a m^cIL cos a 

)+ 0 h 


'n 


'Tit. 


n 


2k +f f 

■] In TTi”™-? -2Ki.cos a coth ^ 


2k -f 
n 




IQI r„ ncj 

] (4.17) 


2kBT 


Now for high energy particles, where k^is large compared to f^, 
we can write In (1 + f/2k^) 4/2k^ and then in the low tempera- 
ture limit coth (hcf/2kgT)2^ one gets after integration, 

2m p'(z.2’,pe^)^ f K, cos a o o o o 

o,yn)= ^^(1+ -V-)[fn- "l-(1-£o/"h] 


me k ^(a^+k, I 


n 


+ 


m^cKi^cos a . f K, cosa o a- /m m o 


m -0-p/T K-^cosa 2 o 

■^e ^ ]_K^k^cosaln(1+ ^)+2[(1 + -^j^) A +2K^k^cos a] 


X [ci(^)cos(^)+si(|^)sin(^)]'|, 


(4.18) 


Mow at the Bragg condition where the two wave picture holds, 
cos a/2k^ = -(Kj^/2k^)^ so that for high energy particles, we 
can neglect cos a/2k^ compared to unity, so that eqn.(4.18) can 

2 


bo written as. 
Re Cj^yn) = 




Mo(A^+Kjj^)B 


f“[l- ^ ln(1+ -^)] 






+ 2fQ^[(!r/ej3)^(1-e-®DA)- e“®:DA]+ ^ K^2^n( l+f^ ^A^) 


2(A"^-K/)[ci(^)cos(^) + si(^)sin(^)]| 


(4.19) 


In the same approximation, eqn. (3.25) for Im C^Q(n) can be 
written as, 
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Txp'(2m )^/^(z. 2 ;^e^)^e ^ f _ , 

Im C^o(n)= ■ — — ^ f^+2f^(-X) ( i_e-®DA) 


Me E.J/^(a2+i£^2^ 


P 


f 

■ A ft; QXi~ 4- n T f A h c \ • /• A ft C \ . c \ ^ Afic \ 

[xan — +ci(^)sin(^)_si(^)cos(^) 


^ksT. 




^kBT- 


•ksT' 


(4.20) 


The Fourier transform of <nl7(r)ln>, the initial state 
expectation value of the interaction potential is given hy, 


^h =(V^‘)/ <3x <nlV(r)|'n> e 
_ V(k^) 


V 


<n|e 


-iK-i,. u 
— n —cr 




n> 


( 4 . 21 ) 


which has hoen calculated hy Kothari and Singwi^^ and one gets, 

2 


^h = 


4tcZ^Z2^ 

■(A^i^) 


p ’ e' 


-W 


(4.22) 


2 


(4.23) 


3h' 


8k. 


B 


mop 

[1+4( ~) (in the low temperature limit) 


The equations (4.19), (4.20) and (4.22) provide the appro- 
priate values of Re Im and 1 -^, to he used in the expres- 
sion (4.11) for intensity. 


4 . 4 D iscussion : 

The extinction distance defined in eqn. (4.12), is 

typically less than 100 i for electron energies above 100 keV 
while the absorption length is quite large (more than 100 

times) compared to This is because Im So is quite small 

compared to (see below). Thus a collection of emitting atoms 
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spread o'vcr several extinction distances from the ciystal sur- 
face.' 5 such as might he achieved hy bombarding the crystal with 
radio-active atoms at a definite energy, will give an intensity 
which is the average of eqn. (4.11) over thickness. This average 
will eliminate the cosine term leaving only the exponentials, 

1 

with the average distance to the emitter replacing t. 

Por typical values of the parameters A = 4.2 x 10 ^ cm“"' 

^n ~ 2 x 1 o”*^ cm"”^ and = lO"*'^, wo see on n-umerical evalua- 
tion that the magnitude of Re is about 10“^ and that of 

^ho^^h about 10~^. Thus the attenuation is very small 
compared to what has been generally chosen phenomenologically . 
In the ideal case, let us study the angular characteristics of 
the intensity pattern for no attenuation (infinite absorption 
length). In this case, from eqn. (4.11), the average intensity 
ie given by, 

<146% = + ( 1 -x)^ 

= 1 - y/(1+y^) (4.24) 

wheie y is approximately given by (from eqn. ( 4 . 12 )), 
sin 20 - 

y (i; — y'T""" (S^-S) ( 4 . 25 ) 

h 

The average intensity <141 > Eas been plotted against y in 
Fig. 1. To include the Bragg reflection on the other side of 
0=0, the curve should be reflected about the value of y, 
corresponding to 0 =0. For electrons, 0=0 corresponds to 
y = -1.5 and for positrons, to y = 1.5. Ibe corresponding 
curves thus obtained have been shown in Figs. 2 and 3. 






t cr^aat 




on 
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It should bo’ noted that the since Re is very small 

compared to V-i^, these curves are essentially identical to those 

36 

oi IcVViiincD ot al. vyithou'c attenuation. Consequently, the 
an^j;iilnr widths of intc'‘nsity patterns arc also equal to the widths 
of their curve s corresponding to no attenuation case. However, 
their cuives with attenuation correspond to an over estimation 
of the attonu'^.tion 0.1) in view of the fact that 

^lio^'^h ^ typical vedues chosen by DeWames et al.^^, 

im])lyin;' thereby very little attenuation. In fact, in our 

calculelion i.s determined onco tho parameters for W^' are fixed 
and need not be chosen phenomenologically. The mass and energy 
dependence of the attenuation (m /B ) shows that for low energy 

^ ir 

hL'avy p'-»rt ides, the attenuation becomes significant. Bor keV 

'7 

electrons and positrons Re 10“^ and the patterns start 

getting, modified duo to this correction. 

From eqn.(4.25) we note- that the thickness averaged inten- 
sity of tho two boom solution has the outer wing in the intensity 
pattern at (y = + 1), 


’W 


+A0g 


■T C' 


■flK 


'h 


(2moBp. 




\hlI^"Sot 2in V2 




( 4 . 26 ) 


1 r ^^h 1 


h p 

(for electrons, at y = 1) 

hi -I "'ho 


V^|-|ReC^,„l 2m^ 1/2 


tlk. 




"h ' -p 

(for positrons, at y = -1) 


(4.27) 
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Thio shows that the condition ^ << 1 for validity of the 

two heniu theory depends upon the mass m^ of the incident particle 

and the strength of the interaction potential V, . 

n 

It also depends upon the particle energy through the 
term Re which is inversely proportional to E^. The mass 
dependence shows that for heavier particles, as soon as A 

jy 

becomes comparable to e^, the two beam theory is not useful. 
Similarly for strong interaction potentials (Tg large) the con- 
dition A e,^/0g <<1 is violated. Actually the particle mass and 
strength of interaction potential, simultaneously govern the 
validity of the two wave picture. Thus, inspite of the large 
mass, (comparable to that of proton) neutrons are describable by 
the two wave picture, since their interaction potential is weak 
compared to that foi proton case (Coulomb or screened Coulomb). 
Similarly, the electrons and positrons, because of their small 
m.ass, may be described by the two beam theory (although marginally, 
because, of strong interaction potential). But for protons, the 
two beam theory is certainly not anolicable. 

The energy dependence appears only when Re ^ho become s 
significant at low energies <?jid under those circumstances, the 
condition A Q-o/Q-a « is more favourably satisfied for positions 
th.an for electrons. Another Interesting consequence of the pre- 
sent theoiy, in this situation, is an indication of the difference 
between widths of electron and positron emission patterns. Erom 
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ein. (4.25) we see that since = (v^ + Oj^^)/Ep, It follows 

that the magnitude of y for e = o Is slightly Inoieased for posi- 
trons and slightly dooreasod for electrons (since Re 0^^ is always 
negatiwo). This shows that whenever Re 0^^ hecomes significant 
(at low energl, s), the angular width of electron emission pattern 
IS mull compared to that of positrons. Although this difference 


IS only qualitative and does not quantitatively explain the obser- 
ved difference^®’ the difference is qualitatively correct and 


in the right direction. This feature is usually absent in the 
two beam theories^® 



CHAPTER V 


ENERGY LOSS OP CffilB&ED PARTICLES 

5.1 Int roduction ; 

In the preceeding Chapters II, III and IV, we have consi- 
dered some of the features of the quantum theory of channeling. 
In particular, we examined the conditions for channeling; energy 
temperature and potential dependence of the phenomena and the 
widths of the emission patterns. One of the most important proh 
lem that needs attention' in the theory of channeling, is the 
problem of energy loss of channeled particles. The quantities 
like the actual penetration depth, the range of the channeled 
particles and the energy spectra of the emitted (or transmitted) 
particles, are essentially determined by the rate of the energy 
loss of the channeled particles. 

Eor neutral particles, such as neutrons, the only mecha- 
nism by which the channeled particles lose energy, is through 
phonon excitations. Eor charged particles, the total contribu- 
tion to energy loss may be divided into three parts; (1) Energy 
loss to phonons (2) Energy loss due to core electrons and 
(3) Energy loss due to conduction electrons. 

The phonon contribution to the rate of energy loss has 
been considered by Lidiard^^. He calculated the rate of energy 
loss due to elastic processes by treating the crystal in 
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hamionic approximaijion. By 2?6garding "the fas1/ particle coming inio 
"bhe ciystial as a pertunbaiionj ii has been shown ihai duning "the 
propagation of particle, phonon shock waves are generated which 
carry away the energy from the particle. For particle velocities 
V, much greater than the velocity of sound, the rat e&f energy 
loss - dk/dt was found to be inversely proportional to V confirm- 
ing the earlier results of Lehman and leibfried^®. Moreover, the 
phonon contribution to energy loss is foimd to be strongly direc- 
tional effect, as expected. 

The energy loss due to core electrons is also expected to 
show directional effects, as was pointed out in Chapter III. The 
energy transferred to core electrons depends upon impact parame- 
ter and must be calculated separately for various core electrons. 
However, a comprehensive theory of this loss has not yet been 
worked out. It is expected that this loss should be small compa- 
red to that due to conduction electrons, because under channeling 

conditions, the particles do not come very close to the ions and 

93 

hence to the various core electrons. In fact Erginsoy-^ has shown 
that it is a good approximation to neglect this and that the energy 
loss to valence electrons is the most important contribution to 
energy loss of charged particles moving under channeling condi- 
tions. 

Thus the stopping power (dE/dx) of charged particles is 
mainly governed by the rate of energy loss due to conduction 
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©lectrons. In 'the first; approximation these constitute a uniformly 
distributed degenerate electron gas and their contribution to the 
total stopping power is not es^ected to show any dependence on the 
direction of motion of the particle in the lattice. This approxi- 
mation is probably better justified in the case of metals, than 
in semiconductors. In the semiconductors, the -valence electron 
density may be nonuniform so that the energy loss does depend 
upon the position of the particle. ■ Actually, this noniiniformity 
of electron density may account for a part of the obser-ved channel- 
ing effects. 


A general formulation of the response of an electron gas 
to a longitudinal disturbance is possible in terms of the complex 
dielectric constant 6(k, m). Such a formulation has been giwen 
by Lindhard^^ and later by lloziores and Pines^^ by using first 
order perturbation theoiy. The essential steps of the treatment 
given by Nozieros and Pines are as follows. 


Consider the unperturbed electron gas, described by the 
Harailttinian, 


li 


P. 


27ce^ 


o 




(5.1) 


i k k^ 

where n is the electron density and is density fluctuations 
given by, 

p, = 2e (5.2) 

JtC 

In the equations (5.1) ani (5.2), summation over i runs over all 
the electrons, whose positions are denoted by x^^. If the Pourier 
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components of the external test charge (incoming particle) are 
written as, e {r^. exp [-i( t + k.r)]+ c.c.} , then the perturbed 
Hamiltonian is given by, 



where n is an infinitesimal (Quantity, How expanding the perturbed 
wave function as, 

-iS Vh 

(t) = S e ^ a (t) 
n 


and calculating the coefficients in first order perturbation theory, 








W + CO 


^'^k'no 


U) -I- 03 

no ^ 


( 5 . 4 ) 


where (*^it^no clement of between the crystal 

eigen states |n> and [ 0> and ^ trans- 
lational invariance, which implies (P ^^^^no ^ 

the perturbed expectation value of P oan be calculated from 
eqn. (5.4), to be (assuming that system is invariant under reflection) 



How using the Pourier expansions of Macroscopic Poisson eqns. 

k.GEj^ = 4716 r^ 6^“"^ 

k.Ej^ = 4TTe[rj^ + <p-^>] 
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One can write, 

i V _ 

. -iort 




I 

‘STFr^jy 


(5.6) 


Then from eqno. (5.5) and (5.6), we get 


Irn 


r(k7aJ)' 


I (5.7) 


Therefore, using the lowest order time dependent perturbation 
theory result for energy loss”'^^ 

,2 2 2 


■»==■#( ^ ) l\\\ K^k^nol ^5.8) 


one can write from oqns. (5.7) and (5.8), 

2 


dW 


871: e I 1 2 _ 1 

^ l^ki SXicTTJJ 


(5.9) 


For high velocities, where Born approximation is satisfactory each 
Fourier component of incoming particle density behaves like a 
test charge of density fluctuation = "I with a frequency u = k.V 
where V is the particle velocity. Therefore the energy loss suffered 
by the fast electron is, 


dE 

TP 



E 

k 



(-•-) lTk7k7vy 


Replacing summation by integration £ = ^ / dk, and performing 

k 

angular integration we get, 

^ - 1^ { ^ ^ ^ er^my” (5.10) 

Thus the final expression for stopping power is given by, 
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dE _ 2^2 r dk ^ , 1 

^ ^ 5 . 11 ). 

If the incident particle has charge z^e, the factor multiplying 
the integration hecomes 2z^^e^/%Y^ . 

Eor high particle velocity, the upper limit of w - inte- 
gration may ho replaced hy and one gets - ^ • But when 

the particle velocity cannot ho chosen as large, the velocity 
dependence of the stopping power is modified depending on the 
foxm of dielectric function. Eor low particle velocities Y << Y^, 
the Eoxmi velocity, the prohlem of stopping power is similar to 
that of the residual resistivity in a dilute metal alloy where 

electrons are scattered hy stationary charged impurities. The 
96 

Mott solution of this prohlem is based on the Bom approximation. 

QA 

According to Lindhard^^ the same perturbation treatment is valid 
for the st opining power problem because the important parameter 
in the collisions is the relative velocity. Eor small Y, the re- 
lative velocity is simply the Eermi velocity Y^,, since only the 
top most electrons can be scattered. 


The first treatment of the stopping power on this basis, 
fox particle velocities small compared to the Eexmi velocity, is 
due to Ecimi and Telle r^'^. They obtained the expression for stop- 


ping power of particle with charge z-jS, as 


4 2-rr 

e m Y 


In ( 


( 5 . 12 ) 


which shows the proportionality of the stopping power to the particle 
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velocity. Hero ra is the electron mass, lindhard^^, lindhard 
axid Winthei and Iruhnikov and Yavlinskii^^ have used the die- 
lectric response ioimulation (5,1) and derived an expression 
which difl'ors only slightly from that of Fermi and Teller, eqn. 
(5.12). For high density electron gas they obtain, 


dE 

Tlx 


o 2 4 2,r 
2z^‘"e^m V 


ixhV. 


[ i3a ( — _ 1 ] 


(5.13) 


Transmission experiments^*^^ in which the energy loss of protons 
and a-particles of 4-30 keV energy has been measured in thin 
films of various matrials, confirm that the stopping power is 
proportional to the particle velocity. 


In the case of heavy ions, nuclear energy losses are 

largei than the inelastic losses to electron gas, upto energy in 

the MeV range. However, nndei channeling conditions, the cross 

section of hat'd nuclear collisions is greatly reduced and it 

b^-coiikis possible to study the electronic stopping power down to 

rel.atively low energies (10-50 keV) where normally inelastic 

losses would be completely masked by nuclear stopping. There 

ii.i increasing evidence that electronic excitation is appreciable 

101 102 

even for low particle velocities. Channeling experiments ’ 
with Xenon ions, over a wide range of energies (1-1000 keV) in 
tungsten single crystals, show that the range of channeled ions 

changes from the R E behaviour, characteristic of low 

^ . 0*5 

velocities, where nuclear stopping power dominates, to the 

behaviour, consistent with an electronic stopping power proporti- 
onal to velocity. A discussion on the velocity dependence of 
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stopping, p.;wer for vnriouL ranges of particle velocity, has been 
recently given by LincLhard”’^^. 


rioxo wo present the results of a calculation loade for 
stopping power using formula (5.11) and different models for 
dielectric function which will give us an insight, as to what extent, 
the stopping power is sensitive to the form of the dielectric func- 
tion. It is seen that for electron velocity upto about ten times 
the Pormi velocity, stopping power is very nearly inversely propor- 
tional to V and is not vury much sensitive to the form of dielectric 
function. It is also seen that for particle velocities small com- 
pared to Pormi velocity, even the use of expression (5.11) leads 
to the result -dE/ciX very nearly proportional to velocity as we see 
from e qn . (5.13). 


5*2 Calculation Using different Approximations for Diele ctric 
Itin ction : 

Por convenience of numerical computation, we write eqn.(5.11) 

in terras of dimensionless variables. Thus expressing k in units of 

2 

k™, the Pebmi raorientum and m in units of bk-/2m, we get, 

o o 


dE 

'33c 


2kK;o 1 

lyr ! ^ s ra Im (-- i ) dm (5.14) 

2 tcK^^ 0^0 e(k,a)) 


with Kq = raV/hkp = V/V^,, being the ratio of particle velocity to the 
Permi velocity and Im [l/€(k,m)] is to bo written in the dimension- 
less units mentioned above. If 6(k, w) = e^(k,w )+ ie 2 (k,“), 

■) 

e^(k,a)) and £ 2 ( 15 :, w) being real, Im e 2 ^^(k,co) 
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ajid since C,^(k,w) vanishes for w > 2k + k^, the integral in 
eqn. (b.14) may be written as, 

2 


2 2(K -1) 

H =, g, r f ° .dk . 
^ 0 0 


wim ( 


-0 


8 ( k , 0) ) 


•) dm 


/ 


2k E 
dk ° 


2(Vi) ’^"o 


J 


Mow the nunierical calculation of the stopping power req.uires the 
Imovvledge of dielectric function. We first give different appro- 
ximatianc in which dielectric function has been calculated. 


(-4) kandorn Phase Approximation (BPA ) : 


The R±^A was developed in the course of studies of screen- 
ing and collective behaviour of an electron liquid. In the early 
work”^^"^, primary emphasis was placed on the explicit introduction 
of collective coordinates (which describe the plasmons) for this 
purpose. The approximation acquired its name on the basis of the 
physical argument that under suitable circumstances, a sum of 
exponential terras with randomly varying phases could be neglected 
compared to N. Thus the approximation, 


P 


k-q 


y p 

Zi 

i 




N 6 


k,q 


was frequently made when 2]£_q a-PP©^'^^®'^ 3- multiplicative factor 

in the basic equations of the theory. This approximation amounts 
to neglecting the exchange diagrams in the G-oldstone diagram 
analysis. Thus the BPl neglects all the exchange correlations 
and represents a good approximation for high density electron 
gas < 1). 
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Tho ( xpression for dielectric function in RPA was first 
calculate d. by lindhard^^, using a self consistent field method. 
This is given by, 


C(k,m) = 1 + Q^(k,o) 


(5.16) 


With, (k,.) = ^ 1, |2i2k±^| 

^ 00 -2k+k^ 

^-pr (5.17) 

2 ~ ^ 

T TC k Thm _ 


(kjco) = 


for oo< 2k-k^ 


'I’T r . /w-k^,2. 


— ^ j for |2k-k^(<oo < 2k+k'^ (5.18) 

4k'^ i-ji 


= 0 for CO > 2k+k'^, 

Q^^^(k,a)) and Q^'^(k,w) representing the real and imaginary parts 
of Q (k,cc). Here k™^ is the Ihomas-I'ermi wave number in units of 


kx, and io given by 


105,106 


3cj 

wMMa|MriA>«rT 

K. p V 




p p 

where co = (47tne /m) is plaana frequency. The fermi wave number 

»P 

is given by”*^^ being defined as AitxJ’/^ = 

0 2 2 
r is generally expressed in Bohr units as r = r a , a = ^ /me 


and for metallic density region r lies between 1.8 and 5.6. Thus 


in EPA wo have , 

i2i ( Trrr 


Qq (k,m) 

~1+QQ^(k,co ) ) ^+( ( k, 


(5.19) 
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(B) Hubbar d Appr o ximation ; 

Hubbard imprcvod upon tbe results of KPA by taking into 
account the shoit- range correlations due to exchange effects. He 
used the dingrt-iniinat ic tcchnicjue to include these effects by suioiiiing 
the exchange diagrams in an approximate way. In this approximation, 
the bxpror.sion for dielectric function is obtained to be, 


C (k,w) 


1 + 


Qo(k, w ) 


1-f(k) Q^(k,m) 


with, 


f(k) = -|- 


r.TT 


1+k"+K 


( 5 . 20 ) 


(5.21) 


s 


whei( ic the screening parameter in units of k-r,. This screen- 
ing psirametur did not appear in the original formula but it was 
added latei by Hubbard himself, in private communication to several 
authors"'^^. Although, usually the screening parameter K_ is re- 

Q 

placed by the Thomas— Fermi wave nimiber k^^, but it has been em- 
phasized"*^®’ that to reproduce the compressibility limit pro- 
perly, K„ should be chosen as given by Noziores-Pines interpola- 

s 

tion formulvU, 

= (4 _ 0.158 k|j,)/(4 + 0.158 k^^,) (5.22) 

2 

For high density metals, this formula yields values of close 
to but somewhat less than k^^^p. Thus in this approximation we get, 

i Qq 

< FT:E;vr > = - [!+{ i-f(k) )Q/(k,.)] -ecc i-f (k) )Qy(i...)] ^ 

(5.23) 
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( C ) K lo inman Aijijroximat ion : 

Ki(;1nmr'n has used self consistent method and 

Intor th. diagrammatic technique’’’''', to calculate the dielectric 

functions t.ilcin/', both exchange and Coulomb correlations into 

accoimt. It has boon shown that Hubbards formula (5.20) has 

correct fonii only in the static limit (os-o) and that even then 

Hubbard’s f(k) is completely incorrect for large k. Actually 

thu static fonnuJa obtained from Kleinman's results which has ' 

boon confiiiuod by variational calculations of Langreth''*^^, is 

found tu give correct results for screening density’'''^ around a 

fixed charged impurity which agree with those obtained by Singwi 
113 

et al. inakin#^: self constent calculations. 


According to Kleinman’s formulation, 




1 


C*(k,-w) 


- 1 - ^ ( 


1 r X (k,id 


^ ' S-TTk; 


X*(k,-a)) 

+ ) (5.24) 


el^^»‘^^ e*^^(k,-w) 


a.x(kj‘’^)+aoX*(lc>-‘’^)+(a.^-a9^)x(k,w) x*(k,-w) 
with ^(k,m) - 1+ — i —A — — 


1 + (a^-a2)x*(k, -w) 


whore a. = 4 


^ — SL1_| ; a _ ^ J 

^ 1+K^’' ^ ^ 1+k‘^+K g 


(5.26) 


and X (k, +w) - x.}(+w) + iX 2 (l“^) where k dependence of the real 
quantities X^(k, +w) and X 2(1^» ±“) been omitted for simpli- 
city, These functions arc given by, 


fr “±^^.2, W+2k±l^~, 

X , ( ) = — -y*' - i f 1 — ( "“5'iy’) ] ~ZZ ? 

^ 2k2 1 Z ' 2k w+2k+k^ 


a)+2k+k^ k^+(D 


k 


(5.27) 


.25) 
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X o(“) 


X o(-w) = 


, 2 

71 k -i-in 


2 2 


ri' r . /w+k-sN 

2i,r I '-(-5F-) ] 

,w_k2,2 


if a3< 2k-k ; 0 otherwise 

(5.28) 


2k 


IJJL I 4 it \ ‘”•1’ . >1 ^ 

5^ j if ic -2k<a)<lc 4'2k; 0 otherwise 


Here C^,^(k,w) given by eqn. (5.25) has significance of being 
diolt-ctric function for an electron. This differs from the die- 
lectric function for a tost charge e(k,u), given by eqn. (5.24), 
since the latter hes no exchange interaction with the responding 
electron gas. xJow using eqns. (5.24) and (5.25), we get. 


Im (1/G(k,m)) = -(1/2) 


{ X2(“)-X2(-“)+(a-|-Q;2^ [x-](-w) X2( 


X^(‘«)X2(-‘^)] y +( I x(-w) j ^-X2(-“) ix(w) I ^2 


/ 


{ 1+a 


* [xi( w) + X^(-w)]+(a^^-a2^) [X-,(w)x^(-w)+X2(“)X2(-“)] + ( afX2(“) 

“X2(“w)]+^«1^-“2^Hxi(~“ ) X2(“)-X-](w)X2(-“)] (5.29 ) 

and Tra( 1/£^^(k,w))= - a^X 2 (“)-a 2 X 2 (-^)+a 2 (a 2 -a^)X^( w) X 2 (-“)+(a^+a 2 ) 


2 2 2 2 

(a^-a2) X2(<^) 1 x(-w) 1 +(2a^ -a^a2-oc2 )x-j (-<JJ ) X2( “) 

’'[x^(w)+X^(~w)]+(a^^-a2''')[x-|(w)X^(-t*j) + X2('^)X2(-“)]} ^ ^ [x:2(<^) 

2 2'^ ^ ^ f \'^'i 2*1 


/ 


{ 1+a 


-X2(~w)l+(ai -a2 ^ [ X,(-w)X2(“)-X^(w)X2(-^^)] > 

(li) Langroth Approximation ; 


(5.30 ) 


Recently langreth has used a variational method to 
calculate the dielectric function. It is shown that Kleinman's theory 
does not use the Hormitian property of interaction potential. 
Consequently his results are correct only in the static case (<*^=0). 
For finite frequencies, langreth gets, 
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G(k,w) = 1 + 




1+f(k,a)) QQ(k,u)) 
*2 


with 


where 


f(k,6)) = 


. _ 1 
'1 " T" 


c^[XL(t^) + X^ ^(-(u)]+2c2Xl(uj)X^ (-U) 




ir 


1+K, 


C2 ~ T 


1 

2 5' 

1+k^+K/ 


(5.31) 


(5.32) 


(5.33) 


s ■ - -s 

and (k, ), which has been written as X^(w) for simplicity, 

is given by, 

f M-f ') 

Xj^(k,w) = _ y -£ 




(5.34) 


where f^ is Permi-fimct ion. langreth did not calculate Xj^Ck^) 
because he did not use it in any numerical calculation. This 
may be calculated in a similar way as lindhard's KPA calculation 
for dielectric function. The detailed steps to calculate the 
real and imaginary parts, ^p'|(^) have been given in 

the Appendix, The result is, 


II 


k-c, r -1 u ^ Ci) +k 2 

(+0)) = _ ^ 1+ ^ [(0+20) In co] _ -^ + ^ [^-(-2k") 1 

~~ 2 2 1 
lnlo3+2k+k2|- ]ln|o>+2k+k^| I for k < 2 (5.35) 

.2 


Z 

471 


■+(o-k 


X^2(“) 


4k‘ 
wkp/ISTtk 


+ TO [■'-( 


0) +k^ • 2. 


2k 


) ]ln 


a)+2k+k^ 


w+2k+k' 


^ I >for k > 2 (5*36) 


for < 2k-k 


TO [l-{ for |2k-lc'^|<«.< 2k+k' 

2 


>!l2(-") 


0 

0 


for w > 2k+k 


(5.37) 


(5.38) 
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Now from eqn. (5.31), we car write, 

*• 

lm(1/£(k,a))) = (pV + pV)ri+Q^(l+F^^)_QV]_( 1 +yRQR_j-lQl) 

^ u 0^0 0^^ 0 ' 

[Q^^+Qj(1+p5]]/[[l+Q^(1+P®)-QjpI]^[Qj(1+pK)+Qy]2] (5.39) 

R. T 

where P and P , the real and imaginary parts of f(k,to) in eqn. 
(5.32), are given by, 

^2 ^L1 ^ ^ (5«40) 

I " 2 

(5.41) 

whore M = ( [Xj^^ ( w)+Xj^., (-w)] Rx^ 2 ^ ^ 

5. 3 Results and Discussion ; 

■Va.rious approximations mentioned above, have been used to 
calculate the stopping power -dE/dx of electrons, using eqn. (5. 15). 
The double integration involved in the expression has been done 
numerically, using Simpson's rule. A satisfactory convergence 
was obtained with 20 as upper limit of k- integration and interval 
0.2 for both, k and w , integrations. The results for electron 
velocities 21 ^ and more upto those corresponding to election 
energy about 1.5 keV, have been given in Tables IV, V, VI, VII 
and VIII, for various densities of conduction electrons (in 
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Table lY : Variation of stopping power (-dE/dx) with 





particle 

energy 


= 2 ). 



^0 = 

V 

Energy Yelocit 
^'e (eY) (lO^cm. 

sec - l ) 

L 

clE/dx) 

in units of 10 

4 

ergs cm 

-1 

EPA Hubbard 
screened 

Hubbard 

unscre- 

ened 

Kleinman 

(elec- 

tron) 

Kle inman 
(Test 
charge ) 

langreth 

2 

50. 12 

4.19 

14*161 

13.676 

13.604 

9.442 

12.989 

14.606 

5 

1 12.77 

7.97 

7.969 

7.747 

7.714 

5.087 

7.395 

8.226 

4 

200.47 

8.39 

5.114 

4.988 

4.970 

5.187 

4.782 

5.269 

5 

513.24 

10.49 

3.581 

3.500 

3.488 

2. 196 

3.365 

3.682 

6 

451.06 

12.59 

2.659 

2.603 

2.595 

1.612 

2.508 

2.751 

7 

613.95 

14.68 

2.060 

2.019 

2.013 

1.238 

1.949 

2.113 

8 

801.89 

16.78 

1.648 

1.616 

1.612 

0.984 

1.563 

1.689 

9 

1014.9 

18.88 

1.351 

1.326 

1.323 

0.802 

1.284 

1.384 

10 

1253.0 

20,98 

1.125 

1.105 

1.102 

0.665 

1.070 

1.151 

11 

1516.1 

23.08 

0.936 

0.920 

0.917 

0.553 

0.891 

0.958 

12 

1804.3 

25.17 

0.787 

0.773 

0.771 

0.464 

0.749 

0.805 

15 

2117.5 

27.27 

0.671 

0.659 

0.657 

0.396 

0.638 

0.686 

14 

2455.8 

29.37 

0.578 

0.568 

0. 566 

0.341 

0.550 

0.592 

15 

2819.1 

31.47 

0.504 

0.495 

0.493 

0.297 

0.479 

0.515 

16 

3207.6 

33.57 

0.443 

0.455 

0.454 

0.261 

0.421 

0.453 

17 

3621 .0 

35.66 

0.392 

0.585 

0.384 

0.231 

0.373 

0.401 

18 

4059.6 

37.76 

0*350 

0.343 

0.343 

0.206 

0.332 

0.358 
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Table V: Variation of stopping power (-cnil/dx) with 
particle energy B , (r = 3). 

P o 


K =-^ 

0 Vp 

Bnorgy 

(eV) 

Velocity 
(loQciri. 
se ( ^ ) 


(-dB/dx) in units., of 10 

_4 

i ergs 

-1 

cm 

EPA 

Hubbard 

screened 

Hubbard 

unscre- 

ened 

Kleinman 
(ele ct- 
ron) 

Kle inman 

(Test 

charge) 

langreth 

2 

22.28 

2.79 

8.249 

8.806 

9.474 

6.895 

9.247 

7.861 

3 

50.12 

4.19 

4.789 

5.029 

5.326 

3.651 

5.178 

4.680 

4 

89. 10 

5.59 

3.115 

3.250 

3.416 

2.269 

3.324 

3.065 

5 

139.72 

6.99 

2. 198 

2.285 

2.391 

1.551 

2.330 

2. 169 

6 

200,47 

8.39 

1*642 

1.702 

1.776 

1.139 

1.732 

1.623 

7 

272.87 

9.79 

1.277 

1.321 

1.376 

0.872 

1.343 

1.264 

8 

356.39 

11.19 

1.025 

1.059 

1.100 

0.691 

1.075 

1.015 

9 

4 5 1 .06 

12.59 

0.843 

0.869 

0.902 

0. 563 

0.882 

0.835 

10 

556.87 

13.99 

0.703 

0.725 

0.751 

0.466 

0.735 

0.697 

11 

556.87 

13.99 

0.585 

0.603 

0.625 

0.387 

0.612 

0.580 

12 

801.89 

16.78 

0.492 

0.507 

0.525 

0.325 

0.514 

0.488 

13 

941.11 

18. 18 

0.419 

0.432 

0.448 

0.277 

0.438 

0.415 

14 

1091.5 

19.58 

0.361 

0.372 

0.386 

0.239 

0.378 

0.358 

15 

1253.0 

20.98 

0.314 

0.324 

0.336 

0.208 

0.329 

0.312 

16 

1425.6 

22.38 

0.277 

0.285 

0.296 

0.183 

0.289 

0.274 

17 

1609.3 

23.78 

0.245 

0.253 

0.262 

0. 162 

0.256 

0.243 

18 

1804.3 

25.17 

0.219 

0.225 

0.233 

0. 145 

0.229 

0.217 
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Table VI: Variation of stopping power (-dE/dx) with particle 
energy E-,(r = 4). 

D 


‘^0= 

~ Energy 
h (eV) 

Vclocit 
( 106 ctn. 
se c" ' ) 

T (- 

(iS/dx) in 

-4 1 

units of 10 «rgs cn 


PUPA 

Hubbard 

screened 

Hubbard 

unscre- 

ened 

Kle ininan 
( elec- 
tron) 

Zle inmar 
(Test 
charge ) 

i 

langreth 

2 

12.53 

2.10 

5.647 

5.984 

6.110 

4.242 

5.849 

5.500 

3 

26.19 

3.15 

3.356 

3.499 

3.554 

2.320 

3.390 

3.361 

4 

50. 12 

4.19 

2,204 

2.284 

2,315 

1.465 

2.214 

2.218 

5 

78.31 

5.24 

1.565 

1.615 

1.655 

1.015 

1.568 

1.576 

6 

112.77 

6.29 

1.172 

1.208 

1.222 

0.748 

1.174 

1,182 

7 

153.49 

7.34 

0.915 

0.941 

0.951 

0.577 

0.915 

0.922 

8 

200.47 

8.39 

0.7^6 

0.756 

0.763 

0.459 

0.736 

0.742 

9 

253.72 

9.44 

0.606 

0.622 

0. 628 

0.575 

0.606 

0.611 

10 

313.24 

10.49 

0.506 

0.519 

0.524 

0.311 

0.506 

0.510 

11 

379.02 

11.54 

0.422 

0.432 

0.436 

0.259 

0.421 

0.425 

12 

451.06 

12.58 

0.354 

0.363 

0.367 

0.218 

0.354 

0.357 

13 

529.37 

13.64 

0.302 

0.309 

0.312 

0.185 

0.302 

0.304 

14 

613.95 

14.68 

0.260 

0.267 

0.269 

0. 160 

0.260 

0.262 

15 

704.79 

15.73 

0.227 

0.232 

0.235 

0.139 

0.227 

0.229 

16 

801.89 

16.78 

0. 199 

0.204 

0. 206 

0.122 

0.199 

0.201 

17 

905.26 

17.83 

0. 177 

0. 181 

0.183 

0. 108 

0. 176 

0. 178 

18 

1014.9 

18.88 

0.157 

0. 161 

0.163 

0.097 

0.157 

0.159 

19 

1130.8 

19.93 

0. 141 

0.145 

0.146 

0.087 

0.141 

0. 142 

20. 

. 1252.9 

20.98 

0.128 

0.131 

0. 132 

0.078 

0.127 

0.129 

21 

1381.4 

22.03 

0. 1 16 

0.119 

0. 120 

0.071 

0.116 

0.117 

22 

1516. 1 

23.08 

0. 105 

0.108 

0.109 

0.065 

0.105 

0. 106 

23 

1657.0 

24.13 

0.096 

0.099 

0.100 

0.059 

0.096 

0.097 

24 

1804.2 

25.17 

0.089 

0.091 

0.092 

0.054 

0.088 

0.089 
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Table YII: Variation of stopping power (~dE/dx) with 
particle energy E (r =5). 

P b 


V 

=_ Energy Volocit 
(oV) (108cm. 
sec-- ' ) 


(-dE/dx) in -units of 10 

-6 -1 

ergs cm 

RPA 

Hubbard Hubbard Kleinman Kleinman 
screened unscre- (elec- (Test langretb 

ened tron) charge) 

2 

8.019 

1.68 

4.229 

4.479 

4.600 

3.147 

4.361 

4.219 

5 

18.04 

2.52 

2.560 

2.664 

2.718 

1.742 

2.562 

2.632 

4 

32.08 

3.36 

1.694 

1.751 

1.781 

1.107 

1.685 

1.745 

5 

50. 12 

4.19 

1.207 

1.244 

1.268 

0.770 

1,198 

1.243 

6 

72. 17 

5.03 

0.907 

0.933 

0.946 

0.570 

0.900 

0.933 

7 

98.23 

5.87 

0.709 

0.728 

0.738 

0.440 

0.704 

0.729 

8 

128.30 

6.71 

0.571 

0.586 

0.593 

0.351 

0.567 

0.587 

9 

162.38 

7.55 

0.471 

0.482 

0.488 

0.287 

0.467 

0.483 

10 

200.47 

8.39 

0.394 

0.403 

0.408 

0.239 

0.391 

0.404 

11 

242.57 

9.23 

0.328 

0.336 

0.340 

0.198 

0.325 

0.336 

12 

288.68 

10.07 

0.276 

0.282 

0.285 

0.167 

0.273 

0.283 

13 

338.80 

10.91 

0.235 

0.240 

0.243 

0.142 

0.233 

0.241 

14 

392.93 

11.75 

0-203 

0.207 

0.210 

0.123 

0.201 

0.208 

15 

451.06 

12.59 

0.176 

0. 181 

0.183 

0. 107 

0,175 

0.181 

16 

513.21 

13.43 

0.155 

0.159 

0. 161 

0.094 

0.154 

0. 159 

17 

579.37 

14.27 

0.138 

0.141 

0.142 

0.083 

0.136 

0. 141 

18 

649.53 

15.11 

0.123 

0.125 

0.127 

0.074 

0. 121 

0.126 

19 

723.71 

15.94 

0. 1 10 

0.112 

0.114 

0.067 

0. 109 

0.113 

20 

801.89 

16.78 

0.099 

0. 102 

0.103 

0.060 

0.098 

0. 102 

21 

'884.08 

17.62 

0.090 

0.092 

0.093 

0.054 

0.089 

0.092 

22 

970.29 

18.46 

0.082 

0.084 

0.085 

0.050 

0.081 

0.084 

23 

1060.5 

19.30 

0.075 

0.077 

0.078 

0.045 

0.074 

0.077 

24 

1154.7 

20. 14 

0.067 

0,071 

0.071 

0.042 

0.068 

0.071 
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tTablc VIII; Variation of stopping power (-dE/dx) 

with particle energy (r =6). 

P s 



Energy 

(ovy 

Velocity 
( loScm. 
so 0 “ ^ ) 


(-dE/dx) 

in units of 10 

4 

ergs cm 

-1 

EPl 

Hubbard 

screened 

Hubbard 

unscre- 

ened 

Kleinman 

(elec- 

tron) 

Kle inman 
(Test 
charge ) 

Langreth 

2 

5.569 

1.39 

5,422 

3.548 

3.625 

2.490 

3.449 

3.103 

3 

12.53 

2.10 

2.091 

2.140 

2.173 

1.589 

2.045 

2.035 

4 

22.28 

2,13 

1.388 

1.414 

1.435 

0.887 

1.352 

1.367 

5 

34.80 

3.49 

0.991 

1.008 

1.020 

0.619 

0.965 

0.980 

6 

50. 12 

4.19 

0.746 

0.757 

0.766 

0.458 

0.726 

0.740 

7 

68.22 

4.89 

0.583 

0.592 

0.598 

0.355 

0.569 

0.579 

8 

89.10 

5.59 

0.470 

0.477 

0.481 

0.283 

0.459 

0.467 

9 

112.77 

6.29 

0.388 

0.393 

0.397 

0.232 

0.379 

0.386 

10 

139. 22 

6 .99 

0.324 

0.329 

0.332 

0.193 

0.317 

0.323 

11 

168.45 

7.69 

0,270 

0.274 

0.276 

0.160 

0.264 

0,269 

12 

200.47 

8,39 

0.227 

0.230 

0.232 

0.135 

0.222 

0.226 

15 

235.28 

9.09 

0.194 

0.196 

0. 198 

0.115 

0.189 

0.193 

14 

272.87 

9.79 

0.167 

0.169 

0. 171 

0.099 

0.163 

0. 166 

15 

313.24 

10.49 

0. 145 

0.147 

0.149 

0.086 

0.142 

0.145 

16 

356.40 

11 . 19 

0.128 

0.129 

0. 151 

0.076 

0.125 

0.127 

17 

402.34 

11.89 

0.1 13 

0.115 

0.116 

0.067 

0.111 

0.113 

18 

451 .06 

12.59 

0.101 

0.102 

0, 103 

0.060 

0.099 

0.100 

19 

502.58 

13.29 

0.091 

0.092 

0.093 

0.054 

o'. 088 

0.090 

20 

556.87 

13.99 

0.082 

0.083 

0.084 

0.049 

0.080 

0,081 

21 

613.95 

14.68 

0.075 

0.075 

0.076 

0.044 

0.072 

0.074 

22 

673.81 

15.38 

0, 068 

0.068 

0.069 

0.040 

0.066 

0.067 

23 

736.46 

16.08 

0.062 

0.063 

0.063 

0.037 

0.060 

0.062 

24 

801.89 

16.78 

0.057 

0.058 

0.058 

0.034 

0.055 

0.057 
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motallic density ravage) i.e. for = 2 to = 6. A comparison 

between the results of different approximations has been made 

fox a particular electron density (r^ = 3) in Pig. 4 where _dE/dx 

has been plotted against V \ Tbe graph shows that results of all 

the approximations for dielectric function show essentially the 

same behaviour i.e, the stopping power to be very nearly inversely 

proportional to the particle velocity in the velocity range of 

2Vp to about 10 Vp. Por higher velocities, the velocity dependence 

tends to become 1/V , as expected, ^In Pig. 4, curve 6 corresponds 

to use of eqn. (5.30) in eqn. (5.14)”]. These results establish 

that the velocity dependence of the energy loss due to conduction 

electrons is rather independent of the model chosen for dielectric 

function. The magnitudes of stopping power also do not differ 

appreciably when different dielectric functions are used. This 

result is to be contrasted with the fact that static impurity 

screening is highly sensitive to the paxticular model used for 

1 1 2— 114 

dielectric function where only static dielectric constant 

(w= 0) is required. 

finally, we have given few results of stopping power of keV 
xenon ions in tuingsten targets in Table IX. The velocity of the 
xenon ions is very small compared to the Fermi velocity. According- 
ly , -dE/dx is found to be very nearly proportional to the particle 
velocity (shown in Pig. 5) as also expected from the formula 
(3.12) or (3.13). 
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Table IX: Variation of stopping power (-dE/dx) of 
155 

Xe in t-ungsten crystal. 


i- Energy Velocity (-(3E/dx) in units of 

S' (kel) (lo^cni. 10"^ ergs cm~‘ 

sec- 'I) KPA Hubbard (screened) 


0.010 

0.467 

0.026 

0.198 

0.226 

0.012 

0.673 

0.031 

0.255 

0.293 

0.014 

0.916 

0.036 

0.288 

0.331 

0.016 

1.196 

0.041 

0.325 

0.375 

0.018 

1.514 

0.047 

0.355 ' 

0.408 

0.020 

1.869 

0.052 

0.365 

0.421 

0.022 

2.262 

0.057 

0.423 

0.488 

0.0 2/( 

2.691 

0.062 

0.463 

0.534 


0.0! 0-02 0-03 0.04 0-05 0.06 0*07 0.08 0.09 O 

(Velocity)"' in Units of 10"® Cm"' sec. 
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CHAPTER VI 


CONCLUSION 


Covtjral interesting conclusions regarding the phenomena 
oi ch.annoling of particles in perfect crystals, emerge from 
the rc'Liults of the inTCstigations presented in this work. The 
fact that the present quantum-mechanical theory confirms some 
of the results expected classically, is in itself interesting 
and is In the usual spirit of quantum mechanics. Thus the 
validity of the condition for anomalous transmission, requires 
that the rrifignitude of the reciprocal lattice vector corres- 
ponding to the plauie (or axis) with which the incident particles 
are being diffracted under the Bragg condition ] K^+k^| =| k^| , 
should be smallest possible. For energetic particle (in which 
case is large compared to 'K^) this means that the incident 
particles should enter into the ciystal very nearly parallel to 
the principal (or low index) crystallographic directions or 
planes. This is what was expected classically in view of the 
fact that low index planes and strings are separated by maximum 
distance between consecutive planes or strings and hence give 
rise to open channels through which particles may propagate with 
minimum attenuation. This aspect is implicit In the Lindhard's 
classical treatment. The effects of lattice vibrations are 
found to decrease the channeling effects and the temperature 
dependence appears as Debye-Waller factor, as one should have expected. 



Then, are several significant results ohta.ined in the pre- 
sent fonn'ilism which charge many of the classical notions, alto- 
gether. It is usually argued that for high energy particles, 
since the deBroglie wavelength of the particle is small compared 
to the lattice spacing, there should not be any significant inter- 
ference effects. On the other hand it is well known that the 
phenomena observed in the electron microscope must be described 
in terms of wave interference: Bragg angles and resonance widths 
dominate the intensity patterns even when the electron wavelength 
is very small compared to the lattice spacing of the target ciys- 
tal. To this end, the present quantum-mechanical treatment shows 
that the existence of anomalous effect has nothing to do with the 
wavelength of the incident particle, but rather is determined by 
two criteria, namely that the interaction potentia,! must be suffi- 
ciently weak compared to the particle energy and that the inter- 
action be localized in the vicinity of the lattice sites. If the 
range of the interaction is small compared to the mean square dis- 
placement, as is expected to be case for neutrons, the attenua- 
tion of both the waves and in 'two beam theory, is 
independent of the potential width and is determined only by the 
mean square displacement. On the other hand, for long range 
interaction potentials, the condition for anomalous transmission 
cannot be satisfied. Comparison between the results of two 
potential models used for charged paarticles, the screened Coulomb 
and the Bom-Mayer, shows that the former is more favourable for 
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channeling. Even for equal range, the penetration depth is 
larger for screened Coulomh potential than for the Born— Mayer 
potential. It is also seen that the channeling condition has 
soMu energy dependent teiras. Although these teims are small at 
high particle cnorgits, they hecorne significant to reduce the 
channeling effect, at low particle energies. Similarly, in 
addition to the usual Dehye-Waller factor, a small temperature 
dependence comes in the channeling-parameter (1-£) which does 
not show any significant effect at low t emperatures. 

It is interesting to note that energy dependence (1/E ) 

Jr 

of thw exponents in the equations (3.17), (3.29) and (3.43), for 

j4j^| . V^hen the channeling parameter (1-6) is not exactly zero 

hut a sma.ll quantity i.o. channeling condition is approximately 

satisfied, the wave 4^ will channel the particles hut there 

will ho some attenuation. The .attenuation upto a given distance 

-1 

traversed will he govk.med hy . Thus the distance at which 

a given fraction of the initially channeled particles is present 
in the channel, is proportional to E^. This conclusion is in 
confirmcLtion with the experimentally established fact that 
the thickness into the ciystal at which one half of the initially 
channeled particles have escaped from the channel is proportional 
to the particle energy E^. 

The quantum theory of the directional effects in the 
charged particle emission in the crystals, presented in the 
Chapter IV, also leads to very important conclusions. The 
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macs and tbu potuntiri dependence of the angles 9^ and A in 
equations (4.26) and (4.27) shows that whenever the p.article 
mass or the strength of potential he cones snail, the interfe- 
rence ,and Bragg diffraction effects modi:fy the emitted intensity 
pattern significantly ( since 4 03/63 becomes small and the 
classical orbital picture is no more applicable.) This com- 
pletely modifies the classical notion of particle trajectory 
being independent of mass for given energy and potential. The 
classical limit is regained only when the particle mass becomes 
very l?i.rgo compared to the electron mass (more than 200 to 300 
times electron mass) because in this limit A 0^ becomes compa- 
rable to 9g and two beam theory loses its validity. Moreover, 
in this formulation, the attenuation of the emitted particle 
emerges automatically as a manifestation of inelastic processes 
in terms of ima.ginary part of renormalization matrix and need 
not be chosen phenomenologically. The actual attenuation is 
found to be small compared to what is usually chosen phenomenolo- 
gically and is proportional to (myEp)^/^^. The temperature de- 
pendence of the attenuation appears as Debye-Waller factor. 
Another intorv sting feature of the present formalism is an indi- 
cation of the difference in the widths of electron and positron 
omission patterns even in the two beam theory. The magnitude of 
this difference is determined by the magnitude of Re which 
is relatively small but has an energy dependence E^" . Thus at 
2. low particle energies, one can see that the positron emission 
pattern has more width as compared to that of electron emission 
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pattt.Tjn, ns has hocn experimentally observed^®’ Although, 
the px'cscnt results do not q^uantitavely explain the experimentally 
nbservecl rlifference, they are certainly in the right direction 
and i^ive a q,ualitative explanation which is usually absent in 
the earlier two bean theories^® 

The calculation of energy loss of charged particles pre- 
sentud in Chapter V shows that the contribution of conduction '* 
electrons to the stopping power is inversely proportiona.1 to 
the particle velocity for velocities greater than tho Poimi ve- 
locity upto about 10 times the Femi velocity. For higher panti- 
do velocities, the dependence goes over to 1/V . From Fig. 4, 
it is also clear th.at the velocity dependence of the stopping 
power is rather insonsitivo to the particular approximation used 
for the dielectric function. Tho actual magnitude is also not 
very sensitive to the form of dielectric function. 

Thus we see that the present q_uant-um mechanical treatment 
of channeling phenomena in crystals leads to many new interesting 
features, not expected classically. It is expected that a more 
detailed analysis using many beam theory may yield few more fine 
structures, with the same basic characteristics obtained here 
in the two beam theory. More realistic models for lattice vib- 
rations may also be examined using formulation of Chapter II. 

The problem of energy loss of chaaineled particles requires a 
little more attention. Firstly an elegant formulation for the 
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conisri'b’u.'tiorj ol corij cloctrons "bo onBrgy loss is very much 
'X'Sir'ihlu because it is expected to shov. directional effects, 
in addition tt; giving a correction to energy loss due to con- 
duction electrons. Secondly, as has been pointed out recently 
that even thC' conduction electron contribution is also expected 
to show s.'uiie diiectional effects and stopping power nay be 
slightly different under channeling conditions than in normal 
conditions. This is significant in semiconductors where the 
density of conduction electrons is not uniform. Therefore, it 
is desirable to improve upon the dielectric foimulation of 
uncTfiy lost! due to conduction electrons and to develop a unified 
foimulation fur the energy loss of charged particles due to con- 
duction and core electrons. 
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APPENDIX 


If we write w , p ard k in energy units (i.e. e = etc.) 


the function X^(k,w) may he written as 

XrCk, w) = Xt.(w)+ ix.p(a)) = - E — 


-^P.^ Vk^ 

O 

CO 


£ “■^p-®p+lc * 


.k'^--2£.k+i n 


^p<'’-Vk) 


— 2 

(2tc)'^ (d -K- 22.k+ir] 


(A.1) 


(A.2) 


Now using the property of Fermi- function f at zero temperature, 
where f p = 1 for p < k^ and 0 otherwise, we can write from (A.2), 
the real and imaginary parts as, 

^F 1 

Xth(w) = - f dp / 42 for k<k^ (A. 5)* 

A% kp_k oo-k -2pkx 

1 rr^"^P ^2. dx 

? U P op J « 

A% 0 -I w -k‘^-2pkx 

^F p 1 , 

+ / p'^dp f p2 Ifor k-<k<2k- (A. 4) 

k-kp x^ 0 ) _k‘^-2pkx ^ 

1 dir 

= - — 5 r I P^dp / s — for k>2k^ (A. 5) 

47c 0 -1 w-k -2pkx ^ 

^F + 1 

and 6(w-k^-2pkx) for w£2idCp-k^ (A. 6 a) 

♦ Hears -x^ = (kp^ - P^ - k^)/2pk. 
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~ 17 o f 6(w-k^-2pkx) 

w-k"^ -1 

for 2kkj,-k^<w<2kkp+k^ (1.61)) 

= 0 for w > 2kkj, + k^ 


Xjj-jC-w) is obtained by replacing w by _ w in the expressions 
(A. 3), (1*4) and (1.5) while "becomes zero because of 

the product of Permi functions ^p( when used in angular 

integration with 6(“+k +2p.jc), gives zero. 


Now performing the angular integrations in eqns. 
(1.6), one gets, 


X 


L1 


. w-k^-2kkp P 

(w) = ?5~>y [ / ^ (w-k -p) la,^p dp 

32Ti:^k^ w+k2_2kkp 

+ 2k^ 2 In p dp] for k<2kp 

a)+k‘^+2kkp 

. w-k^+2kkp 2 

r-75-= [/ (p- oH-k ) In pdp] 

32Tt^k^ w-k^-2kkj, 

for k > 2kj, 


(1.3) to 


(1.7) 


( 1 . 8 ) 


The p-integrations may also be done easily and finally we get, 




for k < 2kp 


(1.9) 
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k. 


■P U-k^ 


^7 


4k 




k-n 1 2 2 

P r / to_k 


for k > 2k™ 


0) 


-k^+2kk. 
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w -k _2kk. 
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0) 

^ ^ irar 


for “ < 2kk™-k'^ 


2 

^ [kp^-(^|-;:)'^] for 2kkp-k^<u)<2kkj.+k2 


TOT" 


P 


(A.10) 
(A. 11) 


If we write these equations in the dimensionless units 


used in the text, we get eqns. (5.35) to (5.38). 



